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DIAGNOSTIC SEQUENCING BY A COMBINATION 
OF SPECIFIC n.RAVAGE AND MASS SP i? fTT^OMETRY 

FEELD HF INVENTION 

5 The present invention is in the field of nucleic acid-based diagnostic 

assays. More particularly, it relates to methods useful for the "diagnostic sequencing" 
of regions of sample nucleic acids for which a prototypic or reference sequence is 
already available (also referred to as 're-sequencing'), or which may be determined 
using the methods described herein. This diagnostic technology is useful in areas that 

10 require such re-sequencing in anipid and reliable way: (i) the identification of the 

various allelic sequences of a certain tegion/gene, (u) the scoring of disease-associated 
mutations, (iii) the detection of somatic variations, (iv) studies in the field of 
molecular evolution, (v) the detemiination of the nucleic acid sequences of 
prokaryotic and eukaryotic genomes; (vO identifying one or more nucleic adds in one 

15 or more biological samples; (vii) and determining the expression profile of genes in a 
biological sample and other areas. 

BACKGRniIND OF INV ENTION 
Complete reference genome sequences for a number of model 
20 organisms as well as humans are currently available or are expected to become 

available in ttie near fiiture. A paraUel chaUenge is to characterize the type and extent 
of variation in the sequences of interest because it underlies tiie heritable differences 
among individuals and populations. In humans, the vast majority of sequence 
variation consists of nucleotide substitutions referred to as single nucleotide 
25 polymorphisms (SNPs). DNA sequencing is the most sensitive method to discover 
polymorphisms [Eng C. and Vijg J. et al. Nature Biotechnol. 15: 422-426 (1997)]. A 
growing panel of such sequence variants, together with powerfijl methods to monitor 
them [Landegren U. et al. Genome Res. 8: 769-776 (1998)], is usefiil in linkage 
studies to identify even tiie most subtle disease susceptibiUty loci [Lander E. and 
30 Schoric N., Science 265: 2037-2048 (1994); Risch N. and Merikangas K., Science 273: 
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1516-1517 (1996)], Also, the identification of all (functional) allelic variants will 
require the re-sequencing of particular regions in a large number of samples 
[Nickerson D. et al. Nature Genet. 19: 233-240 (1998)]. Although a number of 
methods to monitor known SNPs have been developed [Landegren U. et aL, Genome 
5 Res, 8: 769-776 (1998)], re-sequencing is likely to be routmely ^plied to secure 

diagnoses of patients. Indeed, in a significant number of disease-associated genes that 
have been surveyed thus fiir, Uterally hundreds or even thousands of different 
mutations have been identified and catalogued. Consequently, sequence determination 
represents the ultinoate level of resolution and may be the prefenred method to monitor 
10 which mutation or combination of mutations, out of a large number of mutations of 
known clinical relevance, is present. 

It would appear that the field of biomedical genetics will rely heavily 
on sequencing technology. Hence, there is a need for advanced sequencing methods 
that are time- and cost-competitive, and at the same time accurate and robust. Recent 
1 5 developments in this area include improvments to the basic dideoxy chain 
termination sequencing method [Sanger et al Proc. Natl Acad. Sci. USA 74: 
5463-5467 (1977); reviewed by Lipshutz R. and Fodor S. et al. Current Opinion in 
Structural Biology 4: 376-380 (1994)], as well as new approaches that are based on 
entirely new paradigms. Two such novel approaches are sequencing-by-hybridization 
20 (SBH) [Drmanac R. et al. Science 260: 1649-1652 (1 993)] and pyro-sequencing 
[Ronagihi M. et al. Science 281: 363-365 (1998); Ronaghi M. et al. Anal Biochem. 
242: 84-89 (1996)]. While the concepts of these approaches have been experimentally 
validated, their ultimate acceptance and usage may depend on the type of application 
- e.g, de novo sequencing, re-sequencing, and genotyping of known SNPs. 
25 Recently, progress has also been made in the use of mass spectroscopy 

(MS) to analyze nucleic acids [Grain, P.F. and McCloskey, J.A., Current Opinion in 
Biotechnology 9: 25-34 (1998), and references cited therein]. One promising 
development has been the application of MS to the sequence detennination of DNA 
and RNA oligonucleotides [Limbach P., Mass Spectrom. Rev. 15: 297-336 (1996); 
30 Murray K., J, Mass Spectrom. 31 : 1203-1215 (1996)]. MS and more particularly. 
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matrix-assisted laser desoiption/ionization MS (MALDI MS) has the potential of very 
high throughput due to high-speed signal acquisition and automated analysis off solid 
sur&ces. It has hem pointed out that MS, in addition to saving time, measures an 
intrinsic property of ttie molecules, and dierefore yields a significantly more 
5 infomiative signal [KOster H. et al.. Nature Biotecknol, 14: 1 123-1 128 (1996)]. 

Sequence information can be derived directly firom gas-phase 
fragmentation [see for example NordhoflFE. et al, J. Mass Spectrom,, 30:99-1 12 
(1995); Little D. et al, J. A. Chem, Soc, 116: 4893-4897 (1994); Wang B. et al, WO 
98/03684 and WO 98/40520; Bldcker H. et al, EP 0 103 677; Foote S. et al, WO 
1 0 98/5457 1 ]. In contrast, indirect methods measure the mass of fragments obtained by a 
variety of methods in the solution phase, i.e., prior to the generation of gas phase ions. 
In its simplest form, mass analysis replaces the gel-electrophoretic fractionation of 
the fragment-ladder (ic, a nested set of fragments that share one common endpoint) 
generated by the sequencing reactions. The sequencing reactions need not necessarily 
IS be base-q)ecific because the base-calling may also be based on accurate mass 

measurmient of fragments that tenninate at successive positions and that differ from 
one another by one nucleotide residue. The fragment-ladder can be generated by the 
Sanger method [Koster H. et al. Nature Biotechnol, 14: 1 123-1 128 (1996); Reeve 
M.A., Howe R. P., Schwarz T., U.S. 5,849,542; Koster H.. U.S. 5,547,835; Levis R. 
20 and Romano L., U.S. 5,210,412 and U.S. 5,580,733; Chait B. and Beavis R., U.S. 
5,453,247], by base-specific partial RNA digestion [Hahner S. et al. Nucleic Acids 
Res., 25: 1957-1964 (1997); Koster H., WO 98/20166] or by chemical cleavage [Isola 
N. et al. Anal Chem., 71: 2266-2269 (1999); references cited in Limbach P., Mass 
Spectrom. Rev., IS: 297-336(1996)]. An alternative method consists of analyzing 
25 the ladder genmted by exonuclease digestion from either the 3 or 5 *-end [Pieles U. 
et al. Nucleic Acids Res., 21: 3191-3196 (1993); Koster H., U.S. 5,851,765; Engels J. 
et al, WO 98/45700; Tair G. and Patterson D., WO 96/36986; Patterson D., U.S. 
5,869,240]. 

A severe limitation of both the direct and indirect MS methodologies 
30 under the current performance conditions is the poor applicability to chain lengths 
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beyond -30-50 nucleotides. As a consequence, it has been suggested that the 
prospects for MS lie with DNA diagnostic assays, rather than large-scale sequencing 
[Smith L., Nature Biotechnol, 14: 1084-1087 (1996)]. Given the fact that MS 
represents an exquisite means to analyze short nucleotide fragments, the various MS- 
S based processes that have been described for nucleic acid based diagnostic purposes 
generally involve the derivation and analysis of such relatively short firagments [see 
for example Koster H.. WO 96/29431; Koster H. et al, WO 98/20166; Shaler T. et 
al, WO 98/12355; Kamb A., U.S. 5,869,242; Monforte J. et al, WO 97/33000; Foote 
S.e/fl/., WO 98/54571]. 
10 Some of the MS-based assays have been used for the scoring of 

defined mutations or polymorphisms. Other processes derive multiple 
oHgonucleotide fragments and yield a *mass-fingerprint' so as to analyze a larger 
target nucleic acid region for mutations and/or polymorphisms. The latter MS 
analyses are however considerably less infomaative in that they are essentially 
1 5 restricted to the detection of sequence variations. The methods cannot be applied to 
dii^ostic sequencing of nucleic acids, where the term diagnostic sequencing means 
the unequivocal detomination of the presmce, the nature and the position of sequence 
variations. At best, the measurements confirm the base composition of small 
fragments whose masses are determined with sufficient accuracy to reduce the number 
20 of possible compositional isomers. Also, it will be realized that only certain changes 
in composition (as revealed by shifts in the mass spectrum) can be unambiguously 
assigned to a polymorphism or mutation. A match between the spectrum of the 
interrogated sequence and a reference-spectrum obtained from wild-type sequence or 
sequences known to contain a given polymorphism, is assumed to indicate that the 
25 mterrogated nucleic acid region is wild-type or incorporates the previously known 
polymorphisms, thereby disregarding certain other possible interpretations. 

While most methods in the art do yield sequence-related information, 
they do not disclose that a combination of several dififerent mass spectra, obtained 
after complementary digestion reactions, allows for the effective survey of a nucleic 
30 acid region and provides an unambiguous assignment of both known as well as 
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previously imknown sequence variations that occur relative to a reference nucleic acid 
with a known nucleotide sequence. 

In view of the limitations of the methods described above, the art 
would clearly benefit from a new procedure for the diagnostic sequencing of nucleic 
5 acids that would ovocome the shortcomings of the processes discussed above. 

In comparison with conventional sequencing technology, the gel- 
electrophoretic analysis of firagment ladders, the methods of the preset invention are 
more suited for the simultaneous analysis of multiple target sequences. In general, 
each particular sequence or sequence variant is associated with a distinct set of mass 
10 peaks. Consequently, the sequencing reactions according to the methods of the 

present invention lend themselves readily to (i) multiplexing (i.e., the analysis of two 
or more target non-contiguous target regions from a single biological sample), (ii) the 
analysas of hetoozygous samples, as well as (iii) pooling strategies (i.e., the 
simultaneous g^uefncing of the analogous regions derived from two or more different 
1 S biological samples). 

Because of the multiplex capacity, the present methods can be adapted 
as a tool for the genome-wide discovery and scoring of polymorphisms (e.g., SNPs) 
usefid as maikers in genetic linkage studies. The unambiguous 
identification/diagnosing of a number of variant positions is less demanding than full 
20 sequmcing and, consequently, a considerable number of target genomic loci can be 
combined and analyzed at the same time, especially whea their lengths are kept 
relatively small. The number of markers that can be scored in parallel will depend on 
the level of genetic diversity in the species of interest and on the precise method used 
to prepare and analyze the target nucleic acids, but may typically be in the order of a 
25 few tens to up to 100 with current MS capabilities. The addition of multiplexing to 
the higih-precision and high-speed characteristics of MS constitutes a new marker 
technology that enables the large-scale and cost-effective scoring of several (tens of) 
thousands of maricm. Some aspects oftheqiplication of the present methods to 
genome-wide genotyping are described in Example 5. 
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Sequracing reactions according to the methods of the present invention 
yield, in principle, a discrete set of fragments for each mdividual sequence or 
sequence variant whereas conventional sequence ladders stack on top of one another. 
Therefore, such sequraces or sequence variants can be analyzed even when present as 
S a lesser species. This is a useful quaUty for tiie analysis of clinical samples which are 
often genetically heterogeneous because of the presence of both normal and diseased 
cells or in itself (e.^., cancerous tissue, viral quasi-spedes). Additionally, the ability 
to detect mutations at a low ratio of mutant over wild-type allele makes it practicable 
to pool individual biological samples, a strategy which should pennit a more cost- 
1 0 effective search for genomic sequence variations in a population. 

The present invention rests in part on the insight that integration of the 
data obtained in a set of complementary fingerprints produced by an appropriate set of 
complementary cleavage reactions of the invention represents a level of 
characterization of a sample nucleic acid essentially equal to sequence determination. 
1 S The present invention is also directed to the use of cleavs^e protocols that result in the 
generation of cleavage products that range fiom mono- and dinucleotides to fragmoits 
of a few tens of nucleotides that are particularly suited for analysis by MS. At the 
same time, the present method is distinct from the other fragmentation processes that 
are limited to screening target nucleic acids for a wide range of potential mutations. 
20 According to the present invention, a combination of several diffoent mass spectra, 
obtained after conqilementary digestion reactions, coupled with systematic 
computational analysis allows the survey of a selected target nucleic acid or region 
thereof and leads to the unambiguous assignment of botti known and previously 
unknown sequence variations. In certain aspects of the present invention, knowledge 
25 of the reference sequence in combination with the methods disclosed herein allows 
modeling of the experimental sqiproach, anticipation of potential ambiguities, and the 
design of an adequate resolution. 
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SUMMARY OF INVENTION 
The present invention is directed to a mass spectroscopic method for 
detectiiig or analyzing a particular nucleic acid sequence. The present invention is 
useful for de novo sequencing or r&*sequencing nucleic acid in a isqpid and reliable 
5 way which permits, for example, the identification of ttie various allelic sequences of 
a certain region/gene, the identification and scoring of disease-associated mutations, 
the detection of somatic variations, detennining goietic diversity in molecular 
evolution, and the deteraiination of the genomic sequences ey., of viral and bacterial 
isolates. The present invention is also useful for identification of all nucleic acid 
10 molecules in one or more biological samples including for expression profiling i.e., 
identification of all the mRNA species that are expressed in a given cell at a given 
time, by rapidly detennining the sequence of the mRNA that is expressed, 

in one embodiment, the present invention is directed to methods for 
sequence analysis of one or more target nucleic acids for which a known reference 
1 5 nucleic acid sequence is available, in this method, one or more target nucleic acids are 
derived fiom one or more biological sanq>les, and a reference nucleic acid are each 
subjected to complementary cleavage reactions, and the products of tihe cleavage 
reactions are analyzed by mass spectroscopic methods. The mass spectra of the one 
or more target nucleic acids are then compared with the mass spectra of the reference 
20 nucleic acid sequence, and the nucleotide sequence of the one or more target nucleic 
acids is deduced by systematic computational analysis. 

In one aspect of this embodiment, multiple targets, such as cDNA 
clones, are prepared fiom the mRNA of the same biological sample, and are 
separately analyzed as above in parallel experiments. In a second aspect, multiple 
25 targets are derived fiom the same biological sample and are analyzed simultaneously, 
for example in genome-wide genotyping. 

The one or more target nucleic acids may be selected fiom the group 
consisting of a single stranded DNA, a double stranded DNA, a cDNA, a single 
stranded RNA, a double stranded RNA, a DNA/RNA hybrid, and a DNA/RNA 
30 mosaic nucleic acid. 
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In a second embodiment, the one or more target nucleic acids are 
selected from the group consisting of an amplified nucleic acid fragment, a cloned 
nucleic acid fragment, and a series of non-contiguous DNA fragments from the 
genome. In one ssped of this invention, the amplified one or more target nucleic 
5 acids are derived by one or more consecutive amplification procedures selected from 
the groiqp condsting of in vivo cloning, the polymerase chain reaction (PGR), reverse 
transcription followed by fhe polymerase chain reaction (RT-PCR), strand 
displacmient anqilification (SDA), and transcription based piocesses. 

In a preferred embodiment, Ihe amplified nucleic acid fi:agments are 
1 0 RNA transcripts generated from one or more target nucleic acids or a reference 
nucleic acid by a process comprising the steps of: (a) amplifying the one or more 
target nucleic acids or the reference nucleic acid using one or more primers 
cone^nding to a region that is complementary to the one or more target nucleic 
acids or the reference nucleic acid and encoding an e7q)ression control sequence using 
1 5 any one of the amplification procedures described above, and (b) generating RNA 
transcripts fix>m the amplified one or more target nucleic acids or reference nucleic 
acid using one or more RNA polymerases ttiat recognize the transcription control 
sequence on the target or reference nucleic acid The RNA generated by fhe above 
process is then subjected to complementary cleavage reactions to goierate nucleic 
20 acid fragments, which are thm analyzed by MS. The transcription control sequence 
may be selected from the group consisting of an eukaryotic transcription control 
sequence, a prokaryotic transcription control sequence, and a viral transcription 
control sequrace. The prokaryotic transcription control sequence may be selected 
from the groi?) consisting of T3, T7, and SP6 promoters. The cognate RNA 
25 polymerases may be either a wild-type or a mutant form enable of incorporating non- 
canonical substrates with a 2*-substituent other than a hydroxyl groiqi. 

In a third embodiment, die one or more target nucleic acids are 
amplified using modified nucleoside triphosphates. The mass modified nucleoside 
triphosphates may be selected fipom the group consisting of a mass modified 
30 deoxynucleoside triphosphate, a mass modified dideoxynucleoside triphosphate, and a 
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mass modified ribonucleoside triphosphate. The mass modified nucleoside 
triphosphate may be modified on the base, the sugar, and/or the phosphate moiety, and 
are introduced through an eazymaUdo stqp, chemically, or a combination of both. In 
one aspect ibc modification may consist of 2*-substituents other than a hydroxyl group 
S on transcript subunits. In anotfaCT aspect, the modification may consist of 

phosphorothioate intemucleoside linkages or phosphorothioate intemucleoside 
linkages fintfaer reacted wifli an alkylating reagent. In yet another aspect, the 
modification may consist of a methyl group on CS of the uridine-S'-monophosphate 
subunits. Such modifications may alter the specificity of cleavage by certain reagents, 
10 and/or the mass of the cleavage products, and/or the length of the cleavage products. 

In one aspect of the invention, the one or more target nucleic acids and 
reference nucleic acid are subjected to complementary cleavage reactions using 
enzymatic cleavage, chemical cleavage, and/or physical cleavage reactions. In a 
preferred embodiment, fiie one or more target nucleic acids and the reference nucleic 
1 5 acid are subjected to en^matic cleavage reaction using one or more enzymes selected 
fiom the groiq> consisting of endonucleases and exonucleases. In a more preferred 
embodiment, the target nucleic acid is a double-stranded RNA and the endonuclease 
used is a ribonuclease. The ribonuclease may be selected the G-specific T, 
ribonuclease, the A-spedfic ribonuclease, the AAJ specific phyM ribonuclease, the 
20 U/C specific ribonuclease A, the C-specific chicken liver ribonuclease (RNaseCL3), 
and cusativin. In one aspect of this preferred embodiment, the target nucleic acid is a 
phosphorothioate-modified single-stranded DNA or RNA and the endonuclease is 
nuclease PI. 

In another aspect, the mass spectroscopical analysis of the nucleic acid 
25 fragments is performed using a mass spectrometer selected fiY>m the group consisting 
of Matrix-Assisted Laser Desorption/Ionization-Time-of-flight (MALDI-TOF). 
Electro^ray-Ionization (ESI), anid Fourier Transform-Ion Cyclotron Resonance (FT- 
ICR). In a preferred embodiment the mass spectrometer used for the analysis of the 
cleavage fragments is MALDI-TOF 
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In a fifth embodiment, the method of the present invention can be used 
for diagnosing nucleic acid sequence variations in one or more target nucleic acids 
derived from a biological sample, for which a known reference nucleic acid sequence 
is available. In fliis method, one or more target nucleic acids derived from a biological 
S sample, and areference nucldc acid whose sequence has been predetennined are 
subjected to conqilementaiy cleavage reactions, and the products of the cleavage 
reactions are analyzed by mass spectroscopic methods. The mass spectra of the one or 
more target nucleic acids is then compared with the mass spectra of the reference 
nucleic acid, and the nucleotide sequence variations in the one or more target nucleic 
10 acids is then deduced by systematic computational analysis of the sequence variations 
between the one or more target nucleic acids and the reference nucleic acid. A variety 
of acid sequence variations including deletions, substitutions and/or insertions in a 
target nucleic acid can be detramined using the mettiod of the present invention. 

In a sixth enibodimoxt, the method of ttie present invention can be used 
15 for scoring known nucleotide sequence variations in one or more target nucleic acids 
derived from a biological sample, for which a known reference nucleic acid sequence 
is available. In this embodimmt, one or more target nucleic acids derived from a 
biological sample, and a refoence nucleic acid are subjected to complementary 
cleavage reactions, and the products of the cleavage reactions are analyzed by mass 
20 spectroscopic methods. The mass spectra of the one or more target nucleic acid is 
then compared with the mass spectra of the reference nucleic acid sequence, and the 
nucleotide sequence variations^utations in the one or more target nucleic acids are 
scored by comparing the nucleic sequence between the one or more target nucleic acid 
and reference nucleic acid by systonatic computational analysis. 
25 In a seventh embodiment, flie method of the present invention can be 

used for determining the nucleotide sequence {de novo sequencing) of one or more 
target nucldc acids derived from a biological sample for which no reference sequence 
is available. In this method, target nucleic acid, derived from a biological sample is 
subjected to complementary cleavage reactions, and the products of the cleavage 
30 reactions are analyzed by mass spectroscopic methods. The mass spectra of the one 
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or more target nucleic acids coiq[>led with a systematic computational analysis is then 
used to deduce the sequence of the one or more target nucleic acids. 

In an eighth embodiment, the method of flie present invention can be 
used for genome-wide genotyping of one or more known or unknown target nucleic 
5 acids. In this method, one or more target nucleic acids, derived fiom a biological 
sample, are anq^lified and then subjected to complementary cleavage reactions. In one 
aspect, multiple targets are derived fiom a single sample and are analyzed 
simultaneously. The products of the cleavage reactions are then analyzed by mass 
spectroscopic methods. The mass spectra of the one or more known or unknown 
1 0 target nucleic acid is compared with the mass spectra of a reference nucleic acid. This 
comparison is then used to infer the genotype of an organism from which the 
biological sample is derived and to determine therefrom the genetically relevant 
nucleic acid sequence variations of the one or more known or unknown nucleic acids. 

In a ninth embodiment, the mettiod of the present invention can be 
IS used to identify one or more target nucleic acids in one or more biological samples. 
In this method, one or more target nucleic acids, derived from a biological sample, 
are amplified and then subjected to complementary cleavage reactions. In one aspect, 
multiple targets are derived fiom a single sample and are analyzed simultaneously. 
The products of the cleavage reactions are then analyzed by mass spectroscopic 
20 methods. The identity of one or more target nucleic acids is deduced by comparing 
the mass spectra of the one or more known or unknown target nucleic acid with each 
other or by comparison with a plurality of mass spectra of reference nucleic acids. 

In one aspect, the method of the present invention can be used for 
expression profiling, Le. identifying the various mRNA expressed in one or more 
25 biological samples. 

Also encompassed by the present invention is a kit for sequence 
analysis of one or more target nucleic acids using mass spectroscopy, the kit 
comprising a container having one or more sets of reference nucleic acids for which 
the nucleotide sequence is known, one or more nucleic acid cleaving agents, and 
30 computer algorithm/software for comparing the mass spectra of the one or more target 
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nucleic acids with the mass spectra of the reference nucleic acid and deducing 
therefirom the nucleic acid sequence of the one or more target nucleic acids, fa one 
embodiment, the nucleic acid cleaving Bgeat in the kit is a chemical agent, fa an 
alternate embodimmt, the nucleic acid cleavmg agent is an enzyme selected from a 
S groupof enzymes consisting of endonucleases and exonucleases. fa a preferred 

embodiment, the endonuclease is a ribonuclease selected from the group consisting of 
the G-specific Ti ribonuclease, the A-specific U2 ribonuclease, the A/U specific phyM 
ribonuclease, the U/C specific ribonuclease A, the C-specific chicks liver 
ribonuclease (RNaseCL3), and cusativin. 

10 

DESCRIPTION OF DRAWINGS 
Figure 1 A (SEQ ID NO: 1) graphically represents the first 120 
nucleotides of exon 5 of human pS3 as well as the fi:agments that would result from 
cleavage of the (+) and (-) strand transcript after G (RNase-Tl) or A (RNase-U2). The 
1 S dotted and fiiU arrows correspond to the resulting <3-meT and ^4-mer cleavage 

products. The arrows from left to right represent fragments firom the (+) strand, while 
the arrows from right to left represCTt fragments firom the (-) strand. The numbers 
mdicate the neutral molecular masses of die ^A-mcr ribonucleotide firagments. The 
calculation assumes that all fragments contam S'-hydroxyl and 3'-phosphate groups. 
20 Figure IB shows the size distribution of the products that result firom 

base-specific cleavage of a 245 nucleotides long exemplary sequence. 

Figure 2 summarizes the results of the mutational simulation analysis 
of a 200-base-pair segment of the HIV protease gene and shows the percentages of the 
mutational changes that can be detected (hatched bars) and mapped (filled bars). The 
25 results were computed for smgle RNase digests of the (+) and (-) strands with 

respectively RNase-Tl (Tl) and RNase-U2 (U2), sqiarately or combined (T1/U2). All 
refers to the analysis with the four diffoent reactions. 

Figure 3 shows the distributions of the number of diagnostic fragments 
obtamed for the mutational simulation analysis of 1.200 base-pair sequence of HIV 
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v/hen using different lengdi segments of respectively 100, 200, 300, and 600 base- 
pairs. 

Figure 4 sununarizes the results of the mutational simulation analysis 
of 1.200 base-pair sequence of HIV and shows Ae percentages of the single 
S nucleotide substitutions that can be detected O^atched bars) and mapped 

unambiguously (filled bars) as a function of the length of the interrogated segmmts. 

Figure 5 (SEQ ID NO: 2 and SEQ ID NO: 3) is a gr^hic 
representation of the pGEND-Zf(+) derived nucleotide sequences used as a model in 
Examples 2 and 4. The regions corresponding to the PGR primers are underlined. Two 
10 PGR products (1 58 and 1012 base-pairs long) were generated Both amplification 
products encompass the phage T7 promoter site; the transcription initiation site is 
indicated with an arrow. The numbering refers to the respective transcripts (118 and 
972 nucleotides). 

Figure 6 is a graphical representation of the MALDI-TOF mass spectra 

15 oftfaeRNase-A cleavage reactioiisofpGEM3-Z^+) derived traiis(^pts. The 

following transcripts were digested: (A) a regular transcript synthesized with rNTPs, 
(B) a transcript in which UMP residues are replaced by dTMP, (C) a transcript where 
UMP is replaced by dUMP, and (D) one that incorporates dCMP instead of CMP. 
Observed masses are indicated above the peaks that match with predicted digestion 

20 products (see Table II). 

Figure 7A (SEQ ID NO: 4 and SEQ ID NO: 5) is a graphical 
representation of PGR products and transcripts used for diagnostic sequencing of the 
RNase-Tl coding region. Two parallel amplification reactions were performed with 
either the upstream or downstream primer tagged to the T7 promoter. The 

25 amplification products allow the transcription of the (+; upper sequence) or (-; lower 
sequence) strand. The underlined region shows ttie appmded T7 promoter site. An 
arrow indicates the transcription initiation site. 

Figure 7B (SEQ ID NO: 6 through SEQ ID NO: 14)diows the position 
and nature of a number of single, double, and triple mutations in RNase-Tl (reference 

30 denotes the wild-type coding region). 
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Figure 8 is a graphical representation of the MALDI-TOF mass spectra 
obtained for RNase-Tl analysis. Four transcripts were digested with RNase-A: (A) 
dU-incorporating transcript of the (+) strand, (B) dC-transcript of the (+) strand, (C) 
dU-transcript of the (-) strand, (D) dC-transcript of the (-) strand The observed 

S masses of predicted peaks are indicated. Presumed double protonated peaks are 

labeled * with the mass of the parental [M+HT peak indicated between parentheses 
(Figure 8B). One of the peaks in Figure 8D (1207.1+G) is best explained by assuming 
the addition of an extra G-residue at flie transcript 3*-end. Figure 8C only shows the 
900-4800 Da mass range; the digestion product of 1 1 124 Da was not detected. 

10 Figure 9 (panels A, B, and C) is a graphical representation of the 

MALDI-TOF mass spectra of the RNase-A cleavage reaction of a pGEM3-Zf(+) 
derived T7-transcript of 972 nucleotides long. The transcript incorporates dCMP 
mstead of CMP residues. The observed masses of the predicted peaks is indicated. An 
asterisk indicates 2'^'-cyclic phosphate reaction intermediates (see Table V). 

15 

DETAILED DESCRIPTION OF INVENTION 
With current capabilities in mass spectroscq[)y, it is inqiractical to 
sequence nucleic acids greater than -50 bases in loigth. Consequently, an impractical 
and cumbersome number of indqpendent sequencing reactions is necessary to cover 
20 the thousands of bases of a gene or other genetic region of mterest. The methods of 
the present invention described below overcome this limitation. At the same time, the 
present method is distinct from the other fragmentation processes that are limited to 
screening target nucleic acids for a wide range of potential mutations. Indeed, the 
appropriate choice of conqilementary cleavage reactions as described herein allows 

25 the determination of the exact location and nature of a genetic variation. Also, it is 
demonstrated herein that coiiq[>utational protocols are an integral part of the described 
method. The methods and algorithms are required to deduce, on the basis of the 
reference sequence(s), the relation between (i) the spectral changes associated with 
one or more cleavage reactions of a given nature, and (ii) the uniquely defined 

30 sequence variations. 
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Sequencing reactions according to the methods of the present invention 
can be multiplexed, Le, used for the simiiltaneoiis analysis of multiple non-contiguous 
target regions [suprd\. Tha:efore, the methods can be ad^ted as a tool for the 
genome-wide discov^ and/or scoring of polymorphisms (e.g. SNPs) usefid as 
S markers in genetic linkage studies. Indeed, it will be recognized that the unambiguous 
identification/diagnosing of a numbo' of variant positions is less demanding than iiill 
sequencing and that consequently a connderable number of target genomic loci can be 
combined and analyzed in parallel, especially when their Imgths are kept relatively 
small. The niunber of markers that can be scored in parallel will depend on the level 
10 of gaietic diversity in the species of intorest and on the precise method used to 

prepare and analyze the target nucleic acids, but may typically be in the order of a few 
tens to up to 100 or more with current MS capabihties. The addition of multiplexing 
to the high-precision and high-speed diaracteristics of MS constitutes a new marker 
technology that enables the large-scale and cost-effective scoring of several (tens of) 
1 S thousands of markers. Some aspects of the ^plication of the presrat methods to 
genome-wide genotyping are described in Example 5. 

The present invention provides a mass spectroscopy (MS) based 
nucleic acid sequencing method thai overcomes some of the drawbacks inherent in the 
prior art. hi contrast to the previously described methods, the methods of the present 
20 invention do not require the generation of a ladder, i.e an ordered set of nested 
nucleic acid firagments characterized by a common end. Rather, the disclosed 
methods rely on a combination of complementary fragmentation reactions and the 
analytical resolution power of MS to improve mass resolution and mass accuracy. 
The present invention is directed to the use of enzymatic cleavage protocols that result 
25 in the generation of cleavage products that range from mono- and dinucleotides to 
fragments of a few tens of nucleotides that are particularly suited for analysis by MS. 
According to the present invention, a combination of several different mass spectra, 
obtained after complementary digestion reactions, coupled with systematic 
computational analysis allows tiie survey of a selected nucleic acid or region thereof 
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and leads to the unambiguous assignment of both known and previously unknown 
sequence variations. 

The present invention is also directed to methods for the diagnostic 
sequencing (also lefeired to as re-sequencing) of all or part of a sample nucleic acid, 

S i.e. the detemiination of the presence, the nature and the location of the sequence 
variations that occur relative to a related known reference sequ»ce. The sequence 
variations may either be previously identified or hitherto unknown. Diagnostic 
sequencing according to the present invention may focus on particular positions in a 
nucleic acid sequence, e.g. when scoring previously known mutations or 

10 polymorphisms. 

The term "mapping", as used herein, will be understood to include both 
the characterization, Le. detmnination of the nature, and the position of the sequence 
variations. 

The terms 'target DNA", "target sequence*', '*target nucleic acid" and 

IS the like, as used herein, refer to the sequence region which is to be sequenced or re- 
sequenced entirely or in part as well as to the nucleic acid material that is actually 
subjected to one or more complementary cleavage reactions. 

The terms '^ference nucleic acid sequence^, ^'related sequence*', 
"^previously known sequence", and the like, refer to a nucleic acid region, the 

20 sequence of which has previously been detennined which corresponds to the target. 
The reference and target sequences may be found to be identical or may diflfer. The 
reference sequence need not derive from the same species. In many applications, 
several different sequence variants will be available as referoice. The differences 
between a target sequence and its reference sequence may be simple (e.g., single 

25 nucleotide substitutions, deletions and insertions; microsatellite polymorphisms) or 
complex substitution, insertion, and deletion of multiple nucleotides). In certain 
situations, one may not know in advance to what reference sequence, if any, the target 
nucleic acid corresponds. In such situations the interrogated target sequence typically 
corresponds to a portion of a (much) larger refoence sequence and/or to one out of a 

30 plurality of different references. 
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The terms •*unambiguous*', **unique", '^unequivocal", and the like, are 
used to indicate that only a single sequence variation or combination of sequence 
variations can explain (he observed mass spectral changes. 

The trans "complementary (cleavage) reactions", "complementary 
5 cleavages" and the like, as used herein, refer to target nucleic acid digestions 

characterized by vaiying specificity [eg., stringent or relaxed mono- and di-nucleotide 
specificity; digestion with a combination of reagoits; partial cleavage] and/or to 
digestion alternative fonns of the target sequence [eg., the complementary (+) and (-) 
strands; incorporation of modified subunits; analysis of variable portions of the target 
10 sequence]. 

The tenns **transcript" and 'transcription", as used herein, refer to the 
synthesis of a nucleic acid polymer by means of an RNA polymerase. In addition to 
canonical subunits (having a 2'-OH group), a transcript may incorporate non- 
canonical substrates (having any other substituent than a hydroxyl groiq> at the 2*- 
15 position). Canonical and non-canonical substrates may contain additional 
modifications. 

The term "genotyping," as used herein, refers to determining the 
genetic constitution^ which is the particular s^ of alleles inherited by the organism as 
a whole, or the type of allele found at a particular locus of interest. 
20 The term "expression profiling," as used herein, refons to method(s) for 

determining the mRNA expression profile of a givra cell or a population of cells at a 
given time under a given set of conditions. 

Nucleotides are designated as follows. A ribonucleoside triphosphate is 
referred to as NTP or rNTP; N can be A. G, C, U or m^ to denote specific 
25 ribonucleotides. Likewise, deoxynucleoside triphosphate substrates are indicated as 
dNTPs, where N can be A, G, C, T, or U. Througjiout the text, monomeric nucleotide 
subunits are denoted as A, G, C, or T with no particular reference to DNA or RNA. 
When necessary, the nature of the nucleoside monophosphates is clarified by the use 
of more ^ecific abbreviations such as U, m*U, CMP, and UMP to refer to 



wo 00/66771 



PCT/EPOO/03904 



-18- 

ribonucleotides and dC, dU, dCMP, dUMP and dTMP to indicate deoxynucleotides. 
Note that T is not an alternative designation for m^. 



Sequencing via non-ordered sets of specific deavagefiragments 

5 The mediods of the present invention allow the interrogation every 

position in a given target sequence without oeating a fiagment-ladder» i.e. a nested set 
of fragments fliat share one common end^int. The method comprises, in part» 
subjecting one or more target nucleic acids to a set of complementary 
mononucleotide' and/or dinucleotide-specific cleavages, the products of which are 

10 analyzed by mass spectroscopy (MS). A preferred method according to the invention 
includes the specific cleavage of the one or more target nucleic acids at each 
nucleotide by way of two or more sqiarate reactions. The digestion products obtained 
in mononucleotide- and dinucleotide-qiecific cleavage reactions such as those 
described herein range from mononucleotides to fragments of a few tras of 

IS nucleotides and are particularly well suited for analysis by MS. This aspect of the 
invention overcomes flie technical limitation of flie short read lengths encountered 
when analyzing fragment-ladders under the current MS performance. The mass 
^ectra obtained with the methods do not provide a simple readout of the sequence. 
Computational approaches provided herein allow the comparative analysis of the 

20 obtained spectra with those known or predicted for the related reference sequence. 

The abiUty to detect and map sequence variants based on the non- 
ordered set of cleavage fragments according to die present invention resides in part in 
the combiruition of die various complementary site-specific reactions. For example, 
one cleavage scheme useftil in the practice of the preset invention makes use of the 

25 mononudeotide-specific ribonuclease-Tl (RNase-Tl , G-specific) and RNase-U2 (A- 
specific; the limited specificity of this enzyme is recognized and will be dealt with 
below). Both purines (A/G) and pyrimidines (C/T) in a target nucleic acid can be 
examined by cleaving an RNA copy of the two complementary strands of a target 
nucleic acid with both enzymes. MS analysis of the fragments generated by only a 

30 single mononucleotide-specific reaction would detect the presence of most sequence 
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variadons but only a minority of the mutations - in essence those affecting the 
nucleotide that is recognized - would also be localized. Since the methods of the 
present invention examine each of the four bases in a given sequence, each of the 
twelve possible nucleotide substitutions result m the loss of one cleavage site and the 

5 concomitant gain of another cleavage site. This principle is illustrated in Table I for 
the KNase-Tl and RNase-U2 cleavage reactions on the two complementary 
transcripts of a hypothetical target nucleic acid. Transitions affect both the KNase-Tl 
and RNase-U2 cleavage patterns of either the (+) or the (-) strand. As can be seen in 
Table 1 , all transversions change the cleavage pattern of both strands of tfie transcript: 

10 they affect either one of the KNase digests on both strands, or the Tl digest of one 
strand and the U2 digest of the complementary strand. In addition to altering two 
cleavage patterns, each single nucleotide substitution also affects the molecular mass 
of one fragment in each of the remaining two digestion reactions (Table I). In 
conclusion, conq>lCTientary cleavage reactions of the present invention results in a 

1 5 high degree of built-in redundancy. Each nucleotide substitution is potentially 
associated with a maximum of ten differences (data points) with respect to the 
reference spectrum. The loss and gain of a cleavage site are associated with botti the 
disaqppearance and appearance of three peaks; two additional peaks undergo a shift as 
a result of a mass difference. In practice, the 1 Da mass difference between C and 

20 U(T) may result in the loss of a significant amount of information (Table I). More 
particularly, in G- and A-specific cleavage reactions, the C/U transitions may go 
unnoticed while the observed mass difference may not be unambiguously assigned to 
a certain transversion. However, in preferred methods of the present invention 
directed to the analysis of RNA target sequences the method makes use of C and/or U 

25 analogs that exhibit more favorable mass differences, thus allowing the unambiguous 
assigmnent of the mass diffoence to a particular transversion. Example 1 and Table I 
illustrate ttiat S-methyluridine is an cxaaaaple of such a usefiil analog [mHJ; R.I. 
CSiemical, Orange, CA; see also to Hacia J. cr ai. Nucleic Acids Res. 26: 4975-4982 
(1998) for the incorporation of m*UTP during in vitro transcription reactions], 

30 Figure 1 A shows, by way of example, a 120-nucleotide segment of 
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exon 5 of the p53 gene as well as a gr^hical representation of the digestion products 
generated by RNase-Tl and RNase-U2 on an RNA copy of each strand. Figure IB 
displays the size distribution of the base-specific digestion firs^ents derived fit>m 
another exemplary sequence and illustrates that mono-, di- and tri-nucleotides are 
5 considerably more numerous than the la^a digestion products. This distribution is 
expected for mononucleotide specific cleavage reactions that goierate firagments with 
an average length of four nucleotides. Contrary to the size distribution, the number of 
different molecular masses that oUgonucleotides can assume rapidly increases with 
the size of the firagment. Because of the constrained composition of digestion products 
10 {e.g. only one G in the case of RNase-Tl), the number of molecular masses of mono-, 
di- and tri-nucleotides is limited to 1, 3 and 6, respectively. Consequently, mono-, di- 
and tri-nucleotides axe often non-informative in the methods of the present invention 
because their number exceeds (he limited mass space. Figure 1 A illustrates that in 
certain parts of the target sequence one of the cleavage reactions produces many small 
1 5 Augments due to an over-representation of the recognized nucleotide and, 

consequently, yields virtually no information. Howev^ using the method of the 
present invention, this problem is mimmized by the con^>lanentary nature of the four 
reactions which ensures that the fragmmts derived fiom the same region by the othor 
digestions (interrogating under-represented nucleotides) are correspondingly larger. 
20 This indicates a basic attribute of the methods of the present invention. Each of the 
four cleavage reactions yields information about a particular mutational alteration {see 
Table I) and, in general, the redundancy in this information enables the identification 
of the mutation (nature and location) even when part of the information is missing 
from the qpectra as described above. 
25 The methods of the present invoition are therefore largely, yet not 

completely, sequence-independent and permits the re-sequencing of virtually any 
variation. Computer simulations of diagnostic sequencing by the present meOiods, 
more particulariy those involving digestion of RNA copies of each strand with the 
RNases Tl and U2, have shown that for target sequences of up to three hundred base- 
30 pairs --90% or more of all possible single nucleotide substitutions are associated with 
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>4 data points. Fewer than 1% of the substitutions do not result in spectral changes. 
More than 95% of all possible single nucleotide substitutions give rise to unique 
spectral changes and can dierefore be unambiguously identified (see Example 1 and 
Figures 3 and 4). 

5 In sununary, deduction of the sequence according to the methods of the 

present invention is based on the integration of the infonnation that resides in a 
complmentaiy set of ^mass-fingerprints* as well as the previous knowledge about a 
related reference sequmce. The relationship between fliis multitude of data allows 
inferring the presence, nature and position of sequence variations in an unambiguous 

10 way. It is illustrative of the method that the derivation of the sequence is not critically 
dependent on the accuracy, i.e., the absolute values of the mass measurements. It is 
rather the cohort ensemble of mass-shifts and appearances/disappearances of 
cleavage sites that uniquely define the sequence. The computer simulations, described 
henrein, assumed a resolution of S Da or 0.1%, a figure which is well above what can 

15 be achieved with state-of-the-art equipment Also, it should be pointed out that the 
determination of the correct base composition is limited anyway to short fragments, 
even in the case of high-precision measurements [e.^., 5-mers in the case of 
unrestrained sequences and if the measuranent has an accuracy of 0.01% or better, 
Limbach P., Mass Spet^om. Rev. 15: 297-336 (1996)]. Other methods in the art, 

20 which involve the accurate mass detennination to assign the correct base composition 
to one or more fragments, will generally permit the detection of most sequence 
variations but not their unequivocal mapping. In these experiments it is generally 
assumed that a certain experimental observation relates to one particular previously 
known sequence variation, ignoring the fact that alternative sequence variations can 

25 explain the same result. 

The present invention ^compasses several additional embodiments 
and aspects desoibed hereinaft^ and certain oftor embodiments will be readily 
a^arent to one of ordinary sidll in the art. 

30 Target nucleic acid preparation and fragmentation 
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(a) Derivation of target nucleic acid and approaches to cleaving with base-specificity 

Nucleic acid molecules can be isolated fiom a particular biological 
sample using any of a number of procedures, which are well-known in the art, the 
particular isolation procedure chosen being 2q[>pn>priate for the particular biological 
5 sample. To obtain an appropriate quantity of isolated target nucleic acid on which to 
perfonn the methods of the present invention, anq>lification of the target nucleic acid 
may be necessary. Examples of sq^ipropriate amplification procedures for use in the 
invmtion include but are not limited to: cloning [Sambn>ok ^ a/.. Molecular Cloning: 
A Laboratory Manual, Cold Spring Haibor Laboratory Press (1989)], polymerase 
10 chain reaction (PCR) [Newton CR- and Graham A., ?CRy BIOS Publishers (1994)] 
and variations such as RT-PCR [Higuchi aL, Bio/Technology 11: 1026-1030 
(1993)] and allele-specific amplification (ASA), strand displacement amplification 
(SDA) [Terrance Walker G. et aL, Nucleic Acids Res. 22: 2670-77 (1994)], and 
transoiption based processes. 
IS One embodiment of the present invration is directed to methods for 

sequencing (re-sequendng, etc.) Nucleic acid conqprising the digestion of an RNA 
copy of each strand of flie target nucleic acid with the RNases Tl and U2. One of the 
advants^es of the method is fbc use of RNA, which exhibits higher sensitivity and 
better stability in MALDI-MS con^Kured to DNA [Hahner S. el a/.. Nucleic Acids Res, 
20 25: 1957-1964 (1997)], Typically, flie first stage of this aspect of the invention 

involves the amplification of the target nucleic acid by PCR or reverse-transcription 
followed by PCR. (RT-PCR) This can be achieved with a pair of dedicated primers 
that incorporate promote sequences as non-annealing 5'-extensions. In a second 
stage, these promoters are used for the specific transcription of the adjacent sequences 
25 including the target sequences. Preferably, the promoter sequences are small and 
permit the in vitro transcription by a single subunit cognate RNA polymerase such as 
those deriving from bacteriophage T7, T3 and SP6. Preferred for use in this aspect of 
the invention are C and/or U analogs that can be incorporated during transcription and 
ttiat exhibit firvorable mass differences [e.g. m*U; st^ra]. The use of PCR primers that 
30 carry different promoter sequraces permits the generation of an RNA copy of both 
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Strands in two parallel strand-specific transcription reactions. Both strands may also 
be transcribed fix)m the same promoter sequence: this requires two parallel 
amplification reactions with only one promote tagged prima:. Alternatively, the in 
vitro transcripts may also be produced fifom sequences cloned in special purpose 
5 vectors such as the pGEM-type vectors available finom Promega (Madison^ WI) which 
contain appropriate promoters. The tfurd step finrther comprises the treatment of the 
resultant RNA transcripts with one or more complementary mononucleotide-specific 
RNases (e.g RNase-Tl and KNase-U2X such tfiat each desired position in the target 
sequence is interrogated. The final stq> in the process consists of the mass- 
10 spectrometric analysis of the RNA fragments resulting fiom the complementary 

cleavage reactions and the comparison of the spectra obtained with those of the known 
refoence sequence. 

Ahemative schemes to prqpare target nucleic acid obtained fiom a 
biological sample and to subject the target sequence to a set of complementary 
IS mononucleotide-specific cleavage reactions are also wittiin the scope of the invention. 
The target nucleic add can be DNA, cDNA, any type of RNA, DNA/RNA hybrid, or 
of mosaic RNA/DNA composition [depending on the ratio of ribo- and 
deoxyribonucleoside triphosphates (rNTP/dNTP) in the synthesis reaction; Sousa R. 
and Padilla R., EMBOJ. 14: 4609-4621 (1995); Conrad F. et al.. Nucleic Acids Res. 
20 23: 1845-1853 (1995)]. The target sequence may also include modifications tfiat are 
either introduced during or after enzymatic synthesis. 

In general, different forms of each target sequence will be prepared so 
as to be able to perform a complementary set of mono-^ecific cleavage reactions. The 
cleavage reactions may be performed enzymatically and/or chemically. The 
25 mononucleotide-specificity of the digestion reactions may reside in the cleaving agent 
(e.g RNase Tl), in the structure of the target nucleic acid, or in a combination of both. 
For exanq[>le, RNase A (specific for both C- and U-residues) can be made 
monospecific by modifications of tiie substrate sequence that block the ribonucleolytic 
action at C or U residues. RNase A cleavage at U residues can in theory be prevented 
30 by chemical modification [Simoncsits A. ei al. Nature 269: 833-836 (1977)]. The 
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enzymatic incorporation of nucleotide analogs, most notably those modified at the 2'- 
hydioxyl group of the ribose is particularly prefmred in the practice of the invention. 
A variety of such analogs have been demonstrated to be substrates for T7 RNA 
polymerase; e,g, 2*-fluoro, 2*-amino [Aunqp H. et al.. Biochemistry 31: 9636-9641 
5 (1992)], 2"-amethyl [Conrad F. et at.. Nucleic Acids Res. 23: 1845-1853 (1995)], as 
weU as r-deoxy NIPs [SousaR. andPadillaR., EMBOJ. 14: 4609-4621 (1995); 
Conrad F. era/.. Nucleic Acids Res. 23: 1845-1853 (1995)]. The above strategy may 
ako be used to improve fhe specificity of certain RNases such as RNase U2 which is 
said to cleave GpN phosphodiester bonds in esctmsive digests (Brownlee G., in 
10 ^^Laboratory Techniques in Biochemistry and Molecular Biolog/* (Work T.S. and 
Work E.. eds,X North-HoUand, Amsterdam, pp 199-200 (1972)]. Mosaic DNA/RNA 
target sequences that incorporate only one specific rNTP and that can be obtained 
quite efficiently with particular mutant polymerases [Sousa R. and Padilla R.» EMBO 
J. 14: 4609-4621 (1995); Gao G. aL, Proc. Natl Acad. Set. USA 94: 407-41 1 
15 (1997); Bonnin A. et aL. J. MoL Biol. 290: 241-251 (1999)], may aUow mono-specific 
cleavages by alkaline treatment or by digestion with a non-specific RNase such as 
RNase-I [Meador J. et al., Eur. J. Biochem. 187: 549-553 (1990)]. 

Alternative strategies to obtain selective cleavage of target sequences 
make use of phosphorothioate chemistry. DNA and RNA polymers wifii 
20 phosphorothioate intemucleoside linkages in fhe Rp stereo-configuration are readily 
synthesized [see Eckstein F., Ann. Rev. Biochem. 54: 367-402 (1985) and references 
cited therein]. Such phosphorothioate linkages can be specificaUy hydrolyzed 
following alkylation [Gish G. and Eckstein F., Nucleic Acids Symp. Ser, pp 253-256 
(1987); Gish G. and Eckstein F., Science 240: 1520-1522 (1988)]. Mono-nucleotide 
25 specific fragmentation according to this aspect of the invention would require the 
synthesis of targets making use of one particular a-thio nucleotide triphosphate 
substrate. Some nucleases (e.g. nuclease PI) cannot hydiolyze Rp phosphorothioate 
diesters; indirect selective cleavage (at a natural phoqihodiester) may thus be obtained 
with target sequences that incorporate three different aS-dNTPs (or oS-rNTPs). 
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(b) Alternative complementary reactions 

The perfonnance of the present sequencing methods will be understood by those 
skilled in the art to be dependent on the following interrelated factors: (1) the length 
of the region to be sequenced, (2) the resolution of the MS analysis, and (3), to some 

5 extmt, the sequence itself. The longer the region of interest and, consequently, the 
larger the number of digestion products, the more important the resolution becomes. 
Also, die length of the region to be sequenced is directly proportional to the number of 
single micleotide substitutions that cannot be unambiguously mi^ped on the basis of 
the four base-specific fiagmentation patterns only (Example 1 ; Figure 4). Some 

10 sequence motife are intrinsically difloicult to sequence. An example of such a sequence 
is CTAGC,C2C3C4C5GATC (SEQ ID NO: 15), where mutations at C, and Cj cannot 
be discriminated from the same type of mutations at C5 and C4, respectively. Another 
such sequence is GAG,A2G3A4GA, whoe G,->A cannot be discriminated from the 
G3->A mutation; similarly, A2->G and A4->G cannot be distinguished. Finally, the 

1 S four mono-nucleotide q»ecific cleavages may also be insufficient to analyze complex 
sequence variations (see discussion below). Most preferably, therefore, the practicing 
of the presoit invention includes a conqiuter-aided simulation of the re-sequ»cing 
strategy of die intended regioiL Such simulation and analysis will reveal possible 
problematic positions in the sequoice and can be used to assess the usefiilness of 

20 certain additional complementary cleavage reactions as countermeasures to overcome 
such sequencing difficulties. 

One such measure consists of dividing the target region and deriving 
two or more (partially overlqyping) segments (e.g., amplicons) from the sample 
nucleic acid rather than sequencing the target region as a whole. In addition to setting 

25 the length, this allows to exert some control over the composition. This would 

abrogate problems arising when tibe region of interest contains a duplicated segment. 
A second measure consists of carrying out one or more alternative or additional 
reactions involving target fragments that incorporate one or more modified 
nucleotides that exhibit different molecular masses such as is described above. Those 

30 of ^1 in the art will know of ttie existence of a wealth of mass-modified nucleotide 



wo 00/66771 



PCT/EPOOm3904 



-26- 

analogs^ many of which are usefiil and can be reconciled with the enzymatic 
procedures of the present method. The nucleotide analogs will differentially affect the 
masses of many of the digestion products and will therefore yield a significantly 
different spectrum that may reveal the required infonnation. The analogs U and m^ 

S [supra] exemplify this. Simulation studies (which model the present invention) have 
indicated that the use of U resolves certam sequence ambiguities observed with m^ 
(data not shown), while overall the lattor nucleotide analog results in considerably 
fewer sequence ambiguities {see Example 1). 

Another option consists of perfomung one or more additional reactions 

10 on the complementary strand. Compared to, for example, a G-specific cleavage of one 
strand, the C-reaction of the complementary sequence will yield a different set of 
fragments characterized by other mass differences. The effect of including reactions 
on the complementary strand of the target sequence is therefore similar to the use of 
nucleotide analogs. 

IS Still another alternative provided by the present inv^tion and which is 

usefiil in obviating the potmtial problems exemplified above includes using reactions 
with alternative specificities of cleavage. For example, partial base-specific cleavage 
can be achieved by changing the reaction conditions or by use of a specially prepared 
target wherein the cleavable and uncleavable {e.g. 2*-modified; supra) forms of one 

20 particular nucleotide occur randomly. Alternatively, instead of partial base-specific 
cleavages, cme or more q;)ecific digestions diaracterized by a greater stringency can be 
performed (e.g. dinucleotide- or relaxed dinucleotide-specificity; see below). The 
digestion of the target sequence, in double stranded DNA form, with restriction 
enzymes is still another alternative provided by the present invention. Double 

25 digestion (ie. a combination of two base-specific cleavages) of target nucleic acid 
alone or in combination with other digestion methods of the present invention also 
represents an informative alternative within the scope of die present invention. 

Another informative option within the scope of the presmt invention 
involves the analysis of truncated target sequences. More specifically, cleavage of 

30 chain terminated sequences prepared, for exan4>le, by incorporation of a particular 3'- 



wo 00/66771 



PCT/EP00m3904 



-27- 

deoxy nucleotide substrate, will yield spectra that contain additional fragments when 
compared to the spectrum of the full target nucleic acid and will consequoitly provide 
additional information that will, in certain cases, allow a more unambiguous 
indemnification of sequence variation. This approach will be particularly useful for 
5 the charactmzaticni of lengQiy digestion products or regions containing complex 
sequmce variations. 

(c) AUemative complementary reactions: cleavage characierized by a greater than 
mononucleotide specificity 

10 

in still anoAer of its embodiments, the mettiod of the present invention 
also includes nucleolytic processes that are characterized by a dinucleotide- or a 
relaxed dinucleotide-specificity. Such stringency of cleavage will facilitate the 
analysis of longer target sequences because the size distribution of the resultant 

1 5 digestion products is even better suited for analysis by MS than fragments with an 
average length of 4 nucleotides that are generated by mononucleotide-specific 
cleavage. Useful in this aspect of the invention are, for example, restriction 
endonuclease reagents capable of cutting DNA at dinucleotide sequences such as 
those described by Mead D. et al., WO 94/21663 ( PCT/US94/03246). RNases that 

20 preferentially hydrolyze pyrimidine-adenosine (CA and UA) bonds have also been 
identified which are useful in the practice of the presoit invention [£. coli RNase-M, 
Cannistraio V. and KameU D., Eur. J, Biochem. 181: 363-370 (1989); as in an 
endoribonuclease isolated from Saccharomyces cerevisiae, Stevois A. et al„ J. 
Bacterial 164: 57-62 (1985); and as is the Enterobacter sp. C-ribonuclease, described 

25 by Marotta C. et al. Biochemistry 12: 2901-2904 (1973)]. As disclosed and 

exemplified in the present invention, the specificity of these enzymes can, if need be, 
essentially be restricted to CA- or UA-bonds by the use of target nucleic acids that 
incorporate dUMP (or dTMP) on the one hand and dCMP on the other hand. 

Stringent or relaxed dinucleotide-specific cleavage may also be 

30 engineered tfarougih the en^matic and chemical modification of the target nucleic 
acid. By way of non-limiting example, transcripts of the nucleic acid of interest may 



i 
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be synthesized with a mixture of regular and a-thio-substrates and the 
phosphorotbioate intemucleoside linkages may subsequently be modified by 
alkylation using reagents such as an alkyl halide (e.g. iodoacetamide, -iodoethanol) or 
2,3-q>oxy-l-propanoL The phosphotriester bonds formed by such modification are 

5 not expected to be substrates for RNases. Using this procedure, a mono-specific 
RNase, sudi as RNase-Tl, can be made to cleave any tibree, two or one out the four 
possible GpN bonds dqiending on whidi substrates are used in the a-tfiio form for 
target preparation. The rqpertoire of dinucleotid6-q)ecific reagmts usefiil in the 
practice of the present invention may be further expanded by using additional RNases, 

1 0 such as RNase-U2 and RNase-A, In the case of RNase-A, the specificity may be 

restricted to CpN or UpN dinucleotides through the enzymatic incorporation of the 2'- 
modified form of the sqipropriate substrates as described above. For example, to make 
RNase-A specific for CpG dinucleotides, a transcript (target) is prepared using the 
foUowing substrates: aS-dUTP, aS-CTP, aS-ATP, and GTP. Thus, using the 

15 indicated methods described herdn, it is possible to engineer all 16 dinucleotide 

specificities. However, not all dinucleotide-q)ecific reagents described herein would 
be required if the complementary strand of the target nucleic acid is included in the 
analysis. 

The strategy outlined above makes it possible to prevent cleavage 
20 within homopolymor tracts (stretches of A's, G's, C's or T's) by an RNase that is 
made specific (or is made specific as described above) for the repeated nucleotide. 
Indeed, incorporation of a particular aS-NTP, followed by alkylation, will selectively 
prevent cleavage within repeated stretches of that nucleotide, allowing cleavage to 
occur at the 3*-side of the last nucleotide in the rq)eat Simulation studies, similar to 
25 those described in Example 1, have identified diis as a particularly usefiil strategy. 
Sequence analysis by digestion of the two oon^lementaiy strands with RNase-Tl and 
RNase-l]2 yielded a 5- to 10-fold reduction in ttie number of ambiguous mutations 
when aS-GMP and aS-AMP were incorporated in flie respective transcripts. These 
studies also suggest that the selective blockage of cleavage within repeats is 
30 accompanied by a relatively small increase in the average length of the digestion 



wo 00/66771 



PCT/EPOO/03904 



-29- 

pioducts, thereby resulting in considerably less loss of infonnation. 

Those of skill in the art will also readily recognize variations or 
alteniatives in certain aspects ofthe fragmentation methods described herein. Such 
alternatives or variations encompassed by the present invention include but are not 
5 limited to: 

1 . the use of oflier or additional RNases (alone or in combination) having similar 
or aUemative specificities; 



2. the use of mutant or chemically modified RNases with usefiil charactoistics vis- 
10 ^-vis the methods of the present invention [see for example, Loverix S. et aL, 

Nature Struct Biol 5: 365-368 (1998) for an RNase Tl mutant that prefers the 
phosphorothioate analog over the natural phosphodiestCT substrate; see also 
Contreras R, and Ficrs W., FEBSLett. 16: 281-283 (1971) for the production of 
limited digests with a chemically modified RNase]; 

15 

3. ttie use of other nucleotide analogs that exhibit diffi^mt masses and/or 
reactivities, including nucleotides fliat incorporate alternative isotopes; and 



4. alternative specific fragmentation methods^ either chemical [Maxam A. and 
20 Gilbert W., Prva Natl Acad. Sci. USA 74: 560-564 (1977); Richterich P. et al, 

Nudeic Acids Res. 23: 4922-4923 (1995)], or enzymatic. 



Multiplex reactions 

In another embodiment, the methods ofthe present invention are 
25 directed to the simultaneous sequence determination of at least two non-contiguous 
regions in a sample nucleic acid. In contrast to traditional sequencing methods that 
generate a fragment-ladd^ (i.& a nested set of fragments diat share a common 
endpoint), the strategies outlined herein are equally useful for multiplex sequencing. 
Multiplex sequencing, according to the present invention, generally involves the co- 
30 amplification of selected regions of target nucleic acids. This can be achieved by 
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using sets of dedicated primer pairs which flank or are co-tenninal with a target 
nucleic acid to be amplified- Alternatively, the preparation of the multiple target 
nucleic acids comprises the concomitant amplification of restriction firagments derived 
from the sample nucleic acid. Some approaches are illustrated and exemplified in 

S Examples. A q[>ecial case ofmultiplex sequencing consists of the simultaneous 
analysis of the two complementary strands of a double stranded target nucleic acid. 

In yet another embodiment, the methods of the present invention can be 
used for the simultaneous sequence detennination of the corresponding target 
region(s) of at least two biological samples. A sequmce variation in one out of a pool 

10 of analogous target nucleic acids may go unnoticed when analyzing conventional 
sequence ladders by means of gel electrophoresis. With the present methods, a 
sequence variation will, as a rule, yield one or more distinct peaks in the various 
complementary mass spectra. This feature should allow the detection of mutations at 
a significantly lower ratio of mutant to wild-type allele and therefore permit the 

IS analysis of larger pools. The ability to pool riders the present methods usefiil for the 
discovery of sequence variations across particular target regions in a given population. 
For this application, typically S-10 samples may be combined. In case the mutations 
have previously been identified, considerably more samples, several tens, can be 
combined. The characteristics that render the present method usefiil for the analysis 

20 of sample pools make the method also effective for the analysis of heteroTjrgous 
samples (i.e., an equimolar mix of two alleles). 

Mass spectrometric methods 

Mass-spectrometric methods usefiil in the practice of the present 
25 invention include ionization techniques such as matrix assisted laser desorption 

ionization (MALDI) and electrospray (ES). These ion sources can be matched with 
various sqparation/detection formats such as time-of-flight (TOF; using linear or 
reflection configurations), single or multiple quadrupole, Fourier transform ion 
cyclotron resonance (FTICR), ion trqp, or combinations of these as is known in the art 
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of mass spectrometry, pLimbach P., Mass Spectrom. Rev, 15: 297-336 (1996); 
Mamy VL, J. Mass Spectrom.,31: 1203-1215(1996)]. 

Because the present methods generally require the analysis of complex 
oligonucleotide fragment mixtures, the MALDI sq>pioach, mostly resulting in singly 
5 charged molecules, is preferred over ES ^ere significant multiple charging will 
further increase the number of spectral peaks. For the desorption/iomzation process, 
numerous matrixAaser combinations can be employed. 

Sequence determination of simple versus complex variations 
10 In another embodiment, the methods of the present invention are 

directed to the diagnostic sequencing of one or more target nucleic acids that, in 
comparison with a related reference nucleic acid, incorporates a sequence variation 
other than a single nucleotide substitution. Such a sequence variation can involve the 
deletion or insertion of one or more nucleotides as well as the substitution of multiple 
IS nucleotides. 

Similar to singile nucleotide substitutions, the insertion or deletion of a 
single nucleotide represents a simple sequence variation whose analysis using 
methods of the present invention is straightforward. Both of these types of sequence 
variations are associated with a characteristic set of (maximum nine) changes in the 

20 four complCTaentaiy mononucleotide-qpecific fragmentation patterns. It will be 
understood that the methods of the present invention, similar to other sequencing 
methods, may not unambiguously locate the point of insertion or deletion when it 
concerns one nucleotide in a stretch of identical nucleotides. This, however, may be 
taken into consideration when performing a computer assisted analysis of whether the 

25 observed spectra relate in a unique way to a specific sequence variant in accordance 
with the practice of the present invention. 

Analysis of a microsatelhte DNA [also referred to as VNTR (variable 
number tandem repeat) or SSR (simple sequence repeat)] represents a special case 
whose analysis is readily achieved using the mefliods of tfie present invmtion. 

30 Although multiple nucleotides are involved with VNTRs or SSRs, the interpretation 



wo 00/66771 



PCT/EPOO/03904 



-32- 

of the spectral changes on the basis of the known reference sequence is rather simple 
and the polymorphism (an altered nmnber of rq>eat units) may readily be 
characterized. 

The methods of the present invention may also be used to analyze more 

5 complex sequence variations such as diose where multiple nucleotides are affected 
either through insertion, deletion, substitution or a combination tfiereof. The analysis 
of a number of double and triple mutants is described below in Example 3d. Multiple 
substitutions witfiin a target sequence are also expected to be accompanied by a 
characteristic number of spectral changes. This number depends on whether the 

10 substitutions are adjoining or separated, as well as on the intervening sequence in case 
the mutations are separated. Single nucleotide substitutions, isolated by a sequence 
that contains at least one A, G, C, and T, are each associated with 10 spectral 
differences as outlined above, hi general, the analysis of complex sequence variants 
will require (elaborate) conq)utatioiial approaches. One possible algorithm involves 

IS the comparison of the experimentally observed spectra with those generated on the 
basis of all possible sequences in the short region to which the sequence variation is 
confined. Such an algorithm will identify the sequence variant or, in case of 
ambiguities, the diffoent matching sequoaces. This procedure illustrates tiiat the 
present methods may be applied to the die novo sequencing of short regions of a target 

20 sequence. It will be recognized that, in practice, the experimental observations will not 
only set the boundaries but will also define the lengdi of the variant region such that 
the algorithm need not consider insertions or deletions. Additional experimentally 
derived information, such as the absence of a particular nucleotide, can further limit 
the sequence space the algorithm has to explore. In particular ^plications, the 

25 complex sequence variants may be previously known and may thus be part of the set 
of reference sequences. In such cases, the e}q)eiimentally observed spectra may be 
directly correlated to those predicted for the reference sequences. There would 
however still be a need to conq>ute wheOier such correlation is unique. The advantage 
of previous knowledge is that the experimental approach can be ad2q>ted such that the 
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output infoimation indeed relates uniquely to the potentially occuiring complex 
sequoice variations. 

Computer algorithm 

5 The present invention, in part, rests on the insight that computational 

analysis of the spectra obtained in a set of complementary cleavage reactions, and 
comparison of these data with the computationally predicted q>ectral changes fiom the 
known reference sequmce, as illustrated herein, is an important step in the 
unambiguous determination of the presence, the nature and the location of sequence 

10 variations. More specifically, the computational jqpproaches to simulate the 

experiment illustrated herein are necessary to deteraiine whether a unique relation 
exists between the spectra obtained and a particular sequence variation. Accordingly, 
one aspect of the present invention contemplates a method which utilizes a computer 
algorithm or method capable of computing the spectral differences resulting firom one 

15 or more nucleotide differences between the target nucleic acid and the reference 
nucleic acid, the mediod and algoriftm comprising subjecting the reference nucleic 
acid and sequence variants* thereof target nucleic acid having nucleotide 
differraces) to the different base specific cleavages to generate oligonucleotide 
fragments, computing the mass of each oligonucleotide fiagment, generating the mass 

20 spectra of the oligonucleotide firagments firom the reference nucleic acid and the 
sequence variants thereof for each of the base specific cleavage reactions, and 
matching these computationally derived mass spectra with the spectra obtained 
experimentally in the different base specific cleavage reactions. 

In one preferred onbodiment the computer algorithm is designed to 

25 systematically compute the spectra of all possible simple nucleotide variations of the 
reference nucleic acid, including but not limited to all possible single nucleotide 
substitutions, deletions and insertions. Since most of the genetic divosity found in 
living organisms involves single nucleotide variations, most of the experimoitally 
observed sequence variations can be identified with the methods and algorithms of the 

30 present invention, meanmg that one or more matches may be found between the 
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observed spectra and the computationally derived mass spectra. In case a unique 
match is found, the sequence variation in the target nucleic acid is unique. When 
more than one match is found between spectra, the sequence variation cannot be 
established unambiguously. 

5 It will be obvious to the person skilled in the art that diffoent 

sqyproadies may be used for performing the computational analy sis, such as, but not 
limited to, performing the computational analysis on the complete reference sequence, 
or performing a serial computational analysis on segments of the reference sequence 
using, for example, a sliding window. The latter appmsch will enable the 

1 0 identification of different sequence variants occurring in different parts of the 
reference sequence. 

In another embodiment, the methods and computer algorithms of the 
present invention are designed to explore all possible nucleotide sequences in a 
limited segment of the reference sequoice. Such methods and algoriflmis may be used 

15 when the preceding q>proach fidls to give a match, demonstrating that the sequence 
variation does not correspond to a simple nucleotide variation in the reference nucleic 
acid. This may be the case when more than one nucleotide change occurs within a 
short region, such that one or more cleavage products contain multiple nucleotide 
alterations. The region corresponding to these cleavage products can then be explored 

20 further by computing the spectra for all possible sequence permutations and 

determining the matching sequence. It is anticipated that givra sufficient computing 
power, such methods and algorithms may be used for de novo sequencing using mass 
spectral data generated according to the present invention. 

25 Applications of the present methads 

The methods of the present invention are particularly well suited for 
n^idly and accurately re-sequencing nucleic acids finom a variety of biological sources 
including, but not limited to, plants, animals, fungi, bacteria and viruses. Re- 
sequencing implies the detection and mapping of both previously known as well as 

30 unknown sequence variations (e.g. mutations and polymorphisms) relative to a related 
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reference sequence. One of the most notable distinctions with respect to conventional 
gel-electrophoretic analysis of fragment ladders, is that gmerally each particular 
sequence (variation) results in a distinct and characteristic set of (mass) peaks. This 
feature makes the present mettiods effective for die reliable scoring of heterozygous 
S samples, die simultaneous sequraicing of multiple target regions from a single 
biological sample {Le., multiplexing^ as well as &e simultaneous analysis of die 
analogous regions from different samples {Le,y pooling). The use of pools of 
individual samples should permit the cost-effective identification of previously 
unknown sequence variations in a population. Hiis aspect of the invention properties 
10 makes the present methods valuable for clinical and public health studies. Very often 
such studies rely on samples (e.g., saliva, blood, swabs, parafBn-embedded tissue, 
biopsy material) that are cellularly and goietically hetmgeneous and, consequently, 
require assays that can detect mutations at a low ratio of mutant over wild-type allele. 

An additional advantage of the present methodology is that it can be 
1 5 tuned (by reducing the number of complementary cleavage reactions) such that the 
diagnostic sequencing is limited to particular positions in a target nucleic acid, a 
feature useful for the unambiguous scoring of previously identified mutations or 
polymorphisms. The processes described herein can be used, for example, to 
diagnose any of the more than 3000 genetic diseases currently known (e.g., 
20 hemophilias, thalassemias, Ducherme Muscular Dystrophy, Huntington's Disease, 
Alzheimer's Disease and Cystic Fibrosis) or genetic defects yet to be identified. In 
addition, certain DNA sequences may predispose an individual to any of a nmnber of 
diseases or conditions such as diabetes, artherosclerosis, obesity, various autoinunune 
diseases and cancer (e.g., colorectal, breast, ovarian, lung). Depending on the 
25 biological sample, the diagnosis for a genetic disease or genetic predisposition can be 
performed either pre- or post-nalally using the methods of the present invention. Re- 
sequencing of nucleic acids derived from infectious organisms using the methods of 
the present invention may reveal the basis of padiogenicity and may also be useful to 
identify the variation(s) that cause drug-resistance. For example, mutations in the 
30 protease/reverse transcriptase region of the human inununodeficiency virus (HIV) 
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have been implicated in the decreased sensitivity towards the antiviral activity of 
protease and reverse transcriptase (RT) inhibitors. The re-sequendng of the nucleic 
acid encoding these viral domains is therefore of special interest to monitor disease 
progression (see Example 1). Similarly, sequencing, according to the present 

5 invention, niay be usefU to determine the antibiotic-resistance phCT^ 

bacteria [e.g. Mycobacterium tuberculosis; Head S. ei al, Mol Cell Probes 13: 81-87 
(1999); Troesch A. et al., J. Clin. Microbiol. 37: 49-55 (1999)], 

In other embodiments, the present mettiods are directed to the 
identification and classification of target nucleic acids. Analyses according to the present 

1 0 invention characterize nucleic acids at a level essentially equal to sequence determination. 
Thmfore, interrogated unknown sequences may be unambiguously identified by 
conq>arison of the obtained mass spectra with those known or predicted for a plurality of 
reference sequences. In this exercise^ novel sequences that have no matching reference 
database sequence may also be found. The use of the methods for expression profiling 

15 (i.e, flie analysis of cDNA libraries) as well as whol&-genome sequmcmg is exenq)lified 
in Exanq>le 6 and 7, respectively. Other supplications include the detemunation of identity 
or heredity {e.g. , paternity or maternity). 



Kits for practicing the invention 

20 Kits for diagnostic sequencing of one or more target nucleic acids in a 

sample are also provided. In preferred embodimmts, such kits comprise one or more 
reference nucleic acids, various reagents for sequence specific cleavage protocols, and 
computer algorithm(s). Such kits may optionally also contain nucleic acid amplification 
reagents. Additionally, the kits may contain reagents for the preparation of modified 

25 nucleic acids, including but not limited lo modified nucleotide substrates. The kits may 
also contain buffers providing conditions suitable for certain enzymatic or chemical 
reactions. In addition, the kits may contain reagents, such as solid supports, for purposes 
of isolating certain nucleic acids and prqiaring nucleic fiagments for mass spectrometric 
analysis. 

30 The foregoing aspects of the invention are illustrative and should not be 
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constraed to limit the invention as set out in the appended claims. Variations in some 
aspects as well as alternative procedures will be readily recognized by one of ordinary 
skill in the art. 

Example 1 describes modeling the diagnostic sequence analysis of a 1200 
S base-pair region of HIV-1 using methods of tiie present invmtion. 

Example 2 desoibes methods for base-specific cleavage by modifying the 
nucleic acid template to be cleaved. 

Example 3 illustrates the diagnostic sequencing of the RNase-Tl coding 
region according to the methods of the present invention. 
10 Example 4 illustrates the analysis of a -1000 base-pair nucleic acid. 

Example 5 illustrates the use of the present invention for genotyping, 
including multiplex genotyping. 

Example 6 illustrates the use of the present invention for transcription 

profiling. 

IS Example 7 illustrates the use of the present invention for whole genome 

resequencing. 

EXAMPLE 1 

Modeling flic Piagnostic Seancnce Analysis of a 1200 Basc-Dair 

20 Region of HIV-1 

The methods of the present invention have been utilized on a 1 200 
base-pair sequence derived firom human immunodeficiency vims type 1 (HIV-1; 
HXB2 isolate; Genbank accession number K03455; position 2161 to 3360). This 
sequence was used as a model in computer simulations to examine the overall 

2S performance of the method, as well as the occurrence of ambiguities. The selected 
region encompasses the entire protease gene and the first '-270 codons of reverse 
transcriptase [conq>are with Hertogs K. et al., Antimicrob. Agents Chemother, 42: 
269-276 (1998)]. The genotyping/re-sequencing of this domain of clinical isolates of 
HIV is of special interest in order to monitor the emergence of drug resistance- 

30 associated mutations. Single as well as multiple changes have been imphcated in the 
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decreased sensitivity towards the antiviral activity of protease and RT inhibitors 
[Hcrtogs K. et al., Antimicrob. Agents Chemother. 42: 269-276 (1998); Schinazi R. et 
al. Int. Antivir, News 4: 95-107 (1996) and rcferraces cited therein]. 

The principal objective of the computer simulation was to examine the 
5 performance of the re-sequendng method for detecting and mapping SNPs. To this 
end we have performed conqiutational simulation analyses in which we have 
systematically mutated each nucleotide one by one in the 1200 base-pair sequence. 
For each mutation we have calculated the molecular masses of the cleavage products 
that would be generated from a ^ven segment of the sequence in the different four 

1 0 RNase digestion reactions, namely \xpon RNase-Tl and RNase-U2 cleavage of the (+) 
and (-) strands. The comparison of ttiese masses with those of the reference cleavage 
products from the original sequence identifies the masses of the diagnostic fragments 
associated with each mutational change, fragments that either appear or disappear 
as a result of the mutation. The underlying assumption in this analysis was that in 

1 5 order to be measurable, the fragment must have a molecular mass diffoent from those 
of the other cleavage products generated in the same reaction. Furtfamnore, we have 
assumed that the resolution of the mass spec analysis is limited to mass differences 
larger than either 5 Da or 0.1%. In other words fragments whose mass difference with 
other fragments in the same digest is smaller than 5 Da or 0.1% were not scored in the 

20 analysis. The quantitative aspects of a mass spectrum {Le. peak heights) were not 
considered in the present simulation study. For each mutational change we have 
computed the number of fragments that are diagnostic for the presence of the 
mutation. Mutational changes were scored as detectable when there was at least one 
diagnostic fragment (showing a spectral change). In addition, we have examined 

25 whether the mutational changes can also be mapped unambiguously. To this end we 
have compared ttie sets of diagnostic fragments associated with each mutation. 
Mutations that yield unique sets of fragments can be mapped unambiguously, while 
mutations that give the same sets cannot be distinguished from one another. 

In a first simulation analysis we have computed the fraction of SNPs 

30 that may be detected and m^ped using respectively 1, 2 and 4 RNase digestion 
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reactions. To this end we have p^oimed a systematic single nucleotide substitution 
simulation on a 200-bas6-pair segment of the HIV sequence. For each of the four 
different RNase digestion reactions [RNase-Tl and RNa5e-U2 cleavage of the (+) and 
(-) strands] we have calculated the numbor of detectable diagnostic fragments and 

5 have analyzed whefiier these filaments are unique for each mutation. The results 
summarized in Figure 2 show that in each of the single RNase digest reactions a large 
fraction (SS% to 85%) of the mutations are detected. In contrast, only a smaU fraction 
(20% to 30%) of these mutational variations can be mapped unambiguously. The 
principal reason is that many different mutational changes result in the same mass 

10 differences. The fraction of mutations that can be mapped increases to around 60% to 
70% when the data of two RNase digest reactions are combined. The further 
combination of the data from the four different cleavage reactions allows 96% of the 
mutational dianges to be positioned unambiguously and illustrates the advantages of 
the methods of the presmt invention. Close inspection of the sequence ambiguities 

IS reveals that about half of these involve C to U (or conversely A to G) transitions. 

Because the difference in molecular mass between C and U residues is only 1 Da, the 
mass difference in the cleavage products of the strand carrying the pyrimidine base is 
too small to be detectable. Consequently one might expect that these mutational 
changes may become detectable when using m^U instead of U. Computational 

20 simulations using m^ on the same 200 base-pair sequence shows that the fiction of 
mutations that can be mapped unambiguously increases to 98%. Consequently all 
further simulations are based on the use of the analog m^. These results demonstrate 
that the four mononucleotide-specific RNase digests are both necessary and sufiicient 
for re-sequencing of most sequences with a high degree of accuracy. 

25 It will be obvious that the quality of the sequences obtained with the 

methods of the invention will be strongly influenced by the size of the sequence 
segments that are examined. Indeed, the larger the size of the segment, the larger the 
statistical chance that certain relevant diagnostic fragments may coincide with other 
cleavage products generated in the same reaction. We have tfa^efore performed a 

30 systematic single nucleotide substitution simulation analysis on the 1,200 base-pair 
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HTV sequence using different size segments, namely 100, 200, 300 and 600 base- 
pairs. In each simulation a total of 3,600 sii^e mutational substitutions was 
analyzed. For each of fhe four different RNase digest reactions both the number and 
the patterns of flie measurable diagnostic firagments were conq[>uted using the 
S detection limits described above. Figure 3 shows the distribution of the number of 
diagnostic fiagments obtained with the 3,600 mutational changes, in the four different 
analyses. The results clearly indicate that a larger percentage of the single nucleotide 
substitutions is associated with fewer diagnostic spectral changes when using larger 
segments of DNA. 

10 In each simulation we determined both the number of detectable SNPs 

as well as the firaction of SNPs that can be mapped unambiguously. The results of the 
computational simulations summarized in Figure 4 show that almost all the 
mutational changes are detected in the four different analyses. Of the 3,600 SNPs, the 
numbo* that escaped detection were respectively 0, 1, 3 and 9 using 1 00 base-pair, 200 

IS base-pair, 300 base-pair and 600 base-pair segments, respectively. In contrast, the 
fraction of mutational variations that can be mapped unambiguously decreases much 
more when using longer segments. While only 1% of the SNPs are ambiguous when 
analyzmg 100 base-pair segments, that fraction increases to ahnost 10% vdth 600 
base-pair segments. Close inspection of the ambiguities shows that the majority of 

20 these involve nearby (often adjacent) pairs of identical bases where the analysis can 
detemiine the nature of the mutation but fails to identify which of the bases is 
changed. 

In conclusion, the results of the simulations show that the methods of 
the invention are effective for re-sequencing and that even large segments may be 
25 used when oidy a limited numb^ of positions need to be analyzed. Also, it q)pear5 
that in most cases a computer-aided simulation study will be essential in the 
experimental design as well as the data interpretation when using the methods of the 
present invention. Most importantly, the simulations will indicate whether spectral 
changes are unambiguously linked to particular sequence variations. 



30 
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EXAMPLE 2 

Base-Speciflc aeavage by Modification of the Template 
The present example illustrates that the specificity of cleavage by a 
nucleolytic reagent may be furth^ confined througjb ttie modification of the target 
5 template such fliat particular phosphodiester bonds resist cleavage. More particularly, 
it is demonstrated tiiat RNase-A» which normally cleaves at the 3*-side of both C- and 
U-residues, becomes mononucleotide-specific when the target incorporates the 2'- 
deoxy analog of one of these nucleotides. A region of the plasmid vector pGEM3- 
Zf(-»-) (Promega, Madison, WI), encompassing the multi-cloning site as well as the 
10 phage T7 promoter sequences, was used as a model {see Figure 5). 

The first step towards the sequence analysis according to the present 
invention involved the amplification of the 158 base-pair test sequence. The reaction 
was carried out in a total volume of 50 1 using 12.5 pmol each of the forward and 
revorse primer, 200 M of each dNTP, 0.25 1 Taq DNA polymerase (5 U/ 1; Promega, 
1 5 Madison, WI), 1 .5 mM MgCl, and a buffer siqiplied with the enzyme. A&cr an initial 
incubation at 94X for 2 min, 40 cycles of the following temperature program wore 
performed: 94**C for 30 sec, 50**C for 30 sec, and 72**C for 15 sec. The sample was 
kept an additional 15 min at 72^C and then chilled. The PCR reaction product was 
purified (High Pure PCR Product Purification Kit; Roche Diagnostics Belgium, 
20 Brussels, Belgium) and subsequently used for transcription of one q>ecific strand. A 
mutant T7 RNA polymerase (T7 R&DNA™ polymerase; Epicentre, Madison, WI) 
with the ability to incorporate both dNTPs and rNTPs was used in the transcription 
reactions. In addition to a transcription with the regular ribonucleotide substrates, one 
reaction was performed where CTP was replaced by dCTP, while in two more 
25 separate transcriptions dtherdljriP or dTTP replaced XJTP. The transcription 
reactions were run in a 50 1 volume containing: 40 mM Tris-Ac (pH 8.0), 40 mM 
KAc, 8 mM spennidine, 5 mM dithiothreitol, 15 mMMgClj, 1 mM of eachrNTP, 5 
mM of dNTP (in these cases the 2q)propriate NTP was excluded), --40 nM DNA 
template (~2 pmol), and 250 units T7 R&DNA™ polymerase. Incubation was 
30 perfonned at 3TC for 2 hours. After transcription, the fiill-lengfh T7 in vitro 
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tnmsciipts (118 nucleotides) were purified by allowing them to anneal to the 5'- 
biotinylated fonn of the complementary reverse PGR primer (Figure 5) followed by 
capture of the biotinylatcd annealing products onto streptavidin-coated magnetic 
beads. To this end, SO pmol biotinylated reverse primer was added to the transcription 

S reactions. The mixtures were first incubated S min at 70^C and, subsequently, -30 
min at room traiperature. Then, a slight excess of Sera-Mag™ strqitavidin magnetic 
microparticles [Seradyn Inc, hidianapolis, IN; resuspended in SOfil of 2 M NaCl, 20 
mM Tris-HCl (pH 8.0), 2 mM EDTA] was added and ttie resultant mixture incubated 
at room temperature for 30 min with agitation. A magnetic particle collector (MPC; 

1 0 Dynal, Oslo, Norway) was used to collect flie beads, remove the supernatant and, 
subsequently, to wash the beads three times with 100 1 100 mM (NH4)3-citrate. The 
beads were finally resuspended in 3 1 25 mM (NHJj-^jitrate containing 0.5 g bovine 
pancreas RNase-A (SO U/mg; Roche Diagnostics Belgium, Brussels, Belgium) and 
incubated at room tempmture for about 30 min to digest the transcripts to 

1 5 completion. 1 1 of this RNase reaction was removed and added to S 1 matrix solution. 
This 1:1 acetonitrileiH^O matrix solution is saturated widi 3-hydroxypicolinic acid 
(-100 mg^ml), and finther contains 2S mM (NHJj-citrate, (occasionally) 2 pmol/ 1 of 
an oligonucleotide serving as an internal standard, and cation-mchange beads in 
(NH4)^-form (Dowex S0W-X2; Sigma, Samt-Louis, MO) to minimize the presence of 

20 sodium and potassium adducts. After incubating the mixture at room temperature for 
1 5 min, 1 1 was put on the sample plate and allowed to dry. Mass spectra were 
collected using a Reflex III mass spectrometer (Bruker Daltonik GmbH, Bremen, 
Germany). 

The RNase-A cleavage products predicted for each of the four 
25 transcripts are shown in Table U. Note that the mass calculation of the predicted 
firagments assumes a 3 '-phosphate group and not the 2*,3*-cyclic phosphate 
intermediate of the cleavage reaction. Overall, the e3q)erimratally obtamed spectra 
(Figme 6) are in excellent agreement with the predictions. The absence of some of 
the smallest 3-mers (Figure 6A and 6C) may be related to flie mass-gate that was 
30 applied to eliminate the non-informative mono- and di-nucleotide digestion products. 



wo 00/66771 



PCT/EPOO/03904 



-43- 

The predicted 3'-proximal fragment TGTTTC (1830, 1 Da) is only poorly ascertained 
in Figure 6C, ttie spectrum deriving fiom the dU-tianscript This result, along 
with other obsm^ations, suggests ttiat fragments wifli a relatively high dU-content are 
detected with a significantly lower sensitivity using the present MS methodology. 

5 The 2817 Da peak in Figure 6D corresponds to the double protonated form of the 
added oligonucleotide. Some of the expected fragments cannot be resolved because 
they have an identical composition. Also, the digestion products of fhe regular 
transcript that differ by one Da only (e.g. the difference between CMP and UMP; 
Table II) cannot be seen as distinct peaks in Figure 6A. In total, the data convincingly 

10 demonstrate that RNase-A behaves as a C-q)ecific RNase when dTTP or dUTP is 
substituted for UTP, and as a U-specific reagent when dC rather than C is 
incorporated into the substrate transcripts. This high level of nucleobase specificity is 
achieved even und^ the over-digestion conditions used in flie present Example. 

The piotocol described in the present Example is illustrative and certain 

15 modifications and variations will occur to the skilled artisan. The inunobilization of 
the transcripts represents an easy means to prepare die mat^al for MS analysis, 
removal of all other reaction components and exchange of Na^ and counter-ions for 
(NH4)^ (note tihat the subsequent RNase digestion does not require any reagents that 
are 'incompatible' with MS). While other methods, such as chromatography, may be 

20 used to prepare the transcripts or the derived digestion products for MS analysis, the 
present method is favorable in that it is readily amenable to automation and high- 
throughput analysis. In repeat experiments, yielding essentially the same results as 
described herein, the transcripts were digested in water and ~15 nanoliter of these 
digests was directly applied onto a Spectrochip™ (Sequenom Inc., San Diego, CA) for 

25 analysis by MALDI-TOF-MS. 



EXAMPLES 

Diagnostic Sequencing of the RNase-Tl Coding Region 
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The present example illustrates the application of the methods of the 
invention to the re-sequencing of a portion of the RNase-Tl coding region. We 
selected the RNase-Tl coding region because of the availability of a collection of site- 
directed mutants [Steyaert J., Eur. J. Biochem. 247: Ml (1997)] which had 
5 previously been sequenced using the classical dideoxy chain termination method. The 
wild-type and mutant sequraces, used in the present exanqple, are shown in Figure 7. 

a. Analysis of the wild-type RNase-Tl sequence 
The experiments were performed essentially as described in Example 2. 
1 0 First, the selected wild-type RNase-Tl target sequences were amplified by PGR with 
the following primers: 

5'-CCGGATATAAACTTCACGAAGACGG (forward) (SEQIDNO:16) 
5*-GATAGGCCATTCGTAGTAGGGAGAGC (reverse) (SEQIDNO:17) 
The resultant ampUcon was subsequently re-amplified using either a forward or a 
1 S reverse primer that incorporates the T7 promoter site as a S* non-annealing extension 
{see Figure 7A): 

S'-TAATACGACTCACTATAGGGCGACTTCACGAAGACGG (forward) 
(SEQIDNO:18) 

5'-TAATACGACTCACTATAGGGCGAATTCGTAGTAGGGAGAGC (reverse) 

20 (SEQIDNO:19) 

Subsequently, each of the resultant promoter-appended amplicons was used as 
template in two separate transcription reactions. The T7 R&DNA polymerase 
(Epicentre, Madison, WI) was used to prq)are transcripts that incorporate dCMP or 
dUMP instead of respectively CMP and IMP (referred to as the dC- and dU- 

25 transcripts). The transcription reactions were carried out as described in Example 2, 
except that each rNTP was present at 2 mM and incubation was performed overnight 
at 37''C. The four fiill-lcngth T7-transcripts were purified by annealing with a 
biotinylated oligonucleotide that matches with the transcript 3'-end (i.e. the 
biotinylated form of either the forward or the reverse PGR primer used in the first 

30 amplification step) and subsequent capture onto streptavidin microparticles. After 
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extensive washing with (NH4)3-citrate, the transcripts were eluted. The beads were 
resuspended in 3 1 of water and kept at QO^'C for 2 min, immediately followed by 
collection of the beads with the magnet and transfer of the supernatant to a fresh tube. 
Then» the obtained amplified target nucleic acids were digested to completion by the 

5 additionofl lof 100 InM(^ffl4)3<it^ate containing RNase-A. Finally, the 
products were analyzed by MALDI-TOF-MS. 

A graphical rq>resaitation of the spectra is shown in Figure 8A-D. The 
predicted degradation products are listed in Table HI. As witti tfie pGEM3-Zf(+) 
transcripts the obtained spectra are in good agreement with the predictions. A few 

10 peaks that are most likely the result of double protonation were also observed (see 
Figure 8B). The T-reaction on the (-) strand suggests the occurrence of transcripts 
with an extra non-template encoded nucleotide at the 3'-end [MiUigan J. ei al. 
Nucleic Acids Res. 15: 8783-8798 (1987)]. Indeed^ in addition to the expected 3'- 
tenninal fragment, a prominent peak is observed that coincides with the same 

1 S fragment containing an extra G-residue (Figure 8D and Table III). The absence of the 
expected 3*-terminal fragment from the C-reaction on the (+) strand (1 153 Da; Figure 
8A) may be explained by this same phenomenon. In this case, cleavage of the 3'- 
extended transcript would occur and result in the 3'-pho5phoiylated (rattier than the 
3'-OH) form of the predicted fragment, a product which would coincide with another 

20 fragment of the same digestion (1233.7 Da; Table m). 

b. Analysis of selected RNase-Tl single point mutations 
Four single nucleotide substitutions were chosen (mutant #1, #2, #3, 
and #4 in Figure 7B). Each of the mutant sequences was analyzed as described for the 
25 wild-type RNase-Tl coding region (Example 3a). The results are sununarized in 
Table IV. Table IV shows, for each mutation, which 5 fragments of the wild-type 
RNase-Tl reference sequence are affected by the mutation as well as tiie 5 fragments 
that are mutation-specific. It also shows which changes are missing, and consequently 
on how many, out of the ten theoretical data points, the mutation identification is 
30 actually based. Spectral changes are missing because they involve fragments that are 
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too small (<3-mer) or not unique. Also, a few firagments were not experimentally 
observed, e.g., one 3-mer as well as the largest fiagments with a mass of >9,8 Kda. 
Of particular interest are the results conceniing mutation #2. These results indeed best 
illustrate the present invention. In fiiis particular case, all four mono-nucleotide 

5 specific cleavage reactions result in die detection of a mutation, i.e. one will notice 
that die sequence differs fiom the wild-type KNase-Tl coding region. However none 
of these reactions, when Xskea alone^ leads to the unambiguous maiqping of the 
mutation. The C-reaction on the (+) strand results in a new fragment of 1947 Da. Not 
only the single nucleotide mutation #2 can explain the creation of such a 6-mer 

10 [composition = A3G(dU)C]. For example, this is also the case for a double mutation 
that converts the sequence CTACTAC into CAAGTAC {see Figure 7); the TAC peak 
will not be lost because of the presence of a third such 3-mer. The T-reaction on the 
(+) strand results in a q>ectrum where the mass of one fragment has increased by S 6 
Da vAica compared to flie reference spectrum. This suggests the replacement of a dC 

IS by a G. Because the cleavage product contains three dC residues, it is not possible to 
position the substitution. The C-readdon on the (-) strand is at first sight the most 
infonnative; a large reference fragment is afTected by the cleavage. The sequence of 
the fragmmt (GTAGiTT— TGxGATQCSEQ ID NO: 20) is however such timt both the 
G|->C and the G2->C mutation can explain the observed products of 98 14 Da and 

20 1 289 Da [composition = GA(dU)C]. Finally, the T-reaction on the (-) strand is the 
least infonnative and the appearance of a peak of 944 Da [A(dC)U] can be ^plained 
in many different ways. An A(dC)U-fragment is, for example, graerated by 
substitution of the T^-residue for a C in the sequence stretch TAT,TT (see Figure 7). 
In conclusion, mutation #2 exemplifies that in some cases the nature and position of a 

25 sequence variation may only be determined by a combination of at least two different 
complementary cleavage reactions. 



c. Analysis of a mixture of wild-type and mutant RNase-Tl sequences 
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The analyses shown in Table IV can be used to simulate experiments 
where equimolar mixtures of the wild-type RNase-Tl sequence and one of the single 
nucleotide substitutions are examined. In such cases, which mimic heterozygotic 
genotypes, the spectra contain a number of novel fiagments in addition to all those 
S derived finom the (wild-type) reference sequace. The characterization and location of 
the mutation/polymorphism is therefore necessarily based on the novel firagments 
only. Unambiguity requires that the novel fiagmmts are sufBcient to uniquely define 
the mutatioiL Those of skill in the art will realize that zygosity determination is 
straightforward using the present methods because each allele is associated with a 
1 0 distinct set of peaks. 

We performed a number of experiments where on particular single 
nucleotide mutant (eg., mutant #3; Figure 7B) was mixed with wild-type RNase-Tl 
such that the mutant allele was preset at the following fractions: 1 :2. 1 :S, 1 : 1 0, 1 :20, 
1 :50, 1 :100, 1 :200, 1 :500 and 1 :1000, the experiment munics the analyses of pools of 
IS samples diaiacterized by different allele frequencies. First, equivalent quantities of 
the wild-type and mutant target sequences were syntfiesized by PCR amplification 
using conditions where the primers are limiting and completely consumed. After 
mixing the two anq>licons in the desired ratios, the material was re-amplified. Then, 
transcripts of the (-) strand were prepared and digested as described above, except that 
20 transcriptions were performed using all four nucleotide triphosphate substrates in the 
ribo-foim (rNTPs) and that cleavage was carried out with KNase-Tl instead of 
RNase-A, Each of the digestion reactions was measured 5 times. Cleavage with the 
RNase-Tl enzyme generates a polymorphic 15-mer fragment which reads: 
AAAUCAAAACCUUCG(SEQ ID NO: 21), where the underlined residue is changed 
25 to A by mutation #3 {refer to Figure 7A and 7B). The mass of the wild-type and the 
mutant fragment is 4807,91 Da and 4830,95 Da, respectively; the mutation causes a 
shift of 23 Da. We found that there was an excellent linear correlation between the 
allele frequoicies and the relative peak heights (RM),97) and fbat the peak associated 
with the mutant allele could stiU be identified with confidence when it represented 5- 
30 10% of the material. It should be noted that m other experiments flie minimum ratio 
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of mutant over wild-type allele that can be detected might be significantly lower. 
Indeed, in the present example, the reliable detection of the 'mutant peak' was 
somewhat encumbered by the occurrence of an extra peak as evidenced by the control 
spectrum recorded for the wild-type target nucleic acid. This extra peak may possibly 

S be attributed to a low level of Na^-adduct of the wild-type fiagment (22 Da mass 
shift). In all, the latter data indicate that homologous targ^ nucleic acids can be 
pooled and analyzed simultaneously; in addition to reveaUng certain sequence 
variations, the methods of the present invention may permit the allele frequencies to 
be estimated among the pool of biological samples. While diagnostic sequence 

10 detemiination as disclosed herein relies primarily on the appearances and 

disappearances of peaks as well as peak shifts, the present example indicates that 
certain quantitative aspects of a spectrum (e.g., peak height and peak area) can be 
included in the sequence analysis and yield complementaiy valuable information. 

15 d, Anafysis ofRNase-Tl multiple mutants 

The methods of the present invention are not limited to the analysis of 
single nucleotide substitutions. Complex variations can also be sequenced. Table IV 
lists the spectral changes that are predicted to be associated with a number of RNase- 
Tl multiple mutants, more particularly double and triple mutants (mutant #5, #6, #7, 

20 and #8 in Figure 7B). As described above, multiple mutants are associated with a 
characteristic number of spectral changes. In the case of multiple substitutions, with 
no deletions or insertions involved, the number of affected reference fragments is 
always identical to the number of novel fragments. For double mutants the number of 
spectral changes ranges from 12, in case the mutations are adjoining (mutant #5), to a 

25 ffiaximum of 20, in case the mutations are separated by a sequence that contains at 
least one A, G, C, and T. In the latter case, the double mutant is to be treated as two 
concurrent but independent single nucleotide substitutions. Triple mutants are 
associated with a tniniiniim of 14 spectral changes (mutant #7). As with single 
nucleotide substitutions, not all die theoretical spectral changes can or may be 

30 observed and part of the mformation will be lost. In the vast majority of the cases 
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however a systematic computational analysis, based on the obtained spectra and the 
reference nucleic acid sequmce(s), can unambiguously identify and locate the 
sequence variations. 

S EXAMPLE 4 

Mass Spectrometric Analysis of a -^lO OQ Base-Pair Region 
The methods of the invoition are designed to overcome the limitation 
of the short read leng&s encountered with camsA MS-based sequaicmg 
methodologies that involve the analysis of fragment-ladders. One can envision tihat, 

10 depending on the appUcation, target regions of several hundred or even a few 1000 
base-pairs can be analyzed. The present example demonstrates that a large number of 
oligonucleotide fragments can be analyzed simultaneously by the methods of the 
present invention and that, consequently, the detection platform does not impose a 
limit on the methodology. 

15 Following the scheme presented in Example 2, a 1012 base-pair region 

of the plasmid vector pGEM3-Zf(+) (Promega, Madison, WI) was amplified and the 
resultant amplicon, subsequently, used for preparation of a 972 nucleotides long in 
vitro T7 transcript (see Figure 5). The transcript incoiporated dCMP instead of CMP 
such that a U-specific cleavage could be performed by RNase-A. The cleavage 

20 products predicted for this transcript, are listed in Table V. Figure 9 shows the most 
relevant parts of the experimentally obtained spectrum. The primary conclusion from 
the experimental data is that complex mono-nucleotide specific digestion reactions, 
consisting of >200 cleavage products, can be analyzed by mass spectrometry. The 
vast majority of the about 67 predicted distinct peaks are readily identified. Only a 

25 few of the 4-mer fi:agments are not or barely detectable. It also ^pears that in the 
present experiment the assignment of some peaks requires the assumption that (at 
least a portion of) certain digestion products contains a 2%3'-cyclic phosphate instead 
of a 3 '-phosphate group. Such peaks differ fix>m the parent peaks by -1 8 Da. It is well 
known that cyclic phosphates result fiom the tiansesterification cleavage reaction and 

30 that these intermediates get hydrolyzed in a slower second reaction step. 
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EXAMPLES 
Genotyp fa g 

The meOiods of the present invention are also useful for the diagnostic 

S sequencing of multiple non-contiguous le^ons of a sample nucleic acid. This lenders 
the present methods useful for the genome-wide discovery as well as the routine 
scoring of polymorphisms (e.g, SNPs) and mutations at multiple loci in genomic 
DNA. Such multiplex genotyping is conceptually no different than re-sequencing; 
both require that alterations are characterized and positioned unequivocally. Similar to 

10 experiments involving a single target sequence described above, a computer 

simulation can be performed to find out which ones of the observed.spectral changes 
is uniquely hnked to particular genomic alterations. Since multiplex genotyping only 
requires ttie identification/diagnosing of a number of variant positions, it will be 
recognized by those of skill in the art that (i) the complexity (i.e. the combined lengOi) 

IS of the multiple target sequences may be significantly greater than in &e case of full 
re-sequencing, and (ii) a single specific cleavage reaction may often suffice for both 
allele and zygosity identification. Applications which involve the use of two 
sequence-specific cleavages that each positively identify one of the two altemative 
forms of a series of bi-allelic SNPs are also possible using the methods of the present 

20 inventioiL For example, many C to T transitions, the most conmion type of point 

mutations and polymorphisms in human, may be easily scored by a combination of C- 
and T(U)-specific reactions. It is worth mentioning that heterozygous samples 
analyzed using gel-electrophoretic sequencing are often difficult to identify with 
confidence. With the methods described herein, the detection of heterozygosity is 

25 unambiguous because of flie presence of both the wild-type and the mutation specific 
set of mass spectral peaks. 



Multiplex goiotyping will generally involve the co-amplification of 
genomic regions. In the case of previously known SN? genetic markers, co- 
amplification of selected loci can be achieved by using dedicated primer pairs [Wang 
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et aL, Science 250: 1077-1081 (1998)]. Alternatively, a more generic approach can be 
adopted for both the discovery and the subsequent routine scoring of a set of SNPs 
where fhe prqparation of target sequences comprises the concomitant amplification of 
multiple short restriction firagments derived from the sample nucleic acid. This 
S 'random j?ampljng' method may be particularly useful with organisms that have a high 
polymorphism content (e.^., more than 1 SNP m 100 base-pairs). This co- 
amplification can be achieved by ligating to the ends of the restriction firagments 
adaptor sequences that incorporate the target sites for a single PCR primer pair. In 
this approach, the average size of the ampUcons must be small such that the majority 
1 0 incorporates ^1 SNP while, additionally, the total number of the ampUcons must be 
sufficiently small so that their combined length is amoiable to analysis by the present 
methods. These requisites can be met by the appropriate choice of restriction 
enzymes and flie use of methods that permit fhe selective amplification of discrete 
subsets of restriction firagmmts [Vos P. €l id.. Nucleic Acids Res. 23: 4407-4414 
15 (1995); Zabeau M. and Vos P., EP 0534858 (1993); Kikuya Kato, Nucleic Acids Res. 
23: 3685-3690 (1995)] and as described herein. For example, a first restriction 
enzyme that cleaves rarely in the genome under study can be combined with a second 
reagent that generates fragments with an average size of about 100 base-pairs (e.g., a 
combination of two en^mes with tetra-nucleotide recognition sites). The number of 
20 fiagments edged by the two different restriction sites should preferably be less than 
100,000; a suitable subset of these can readily be amplified by the use of selective 
primers [Vos P. et aL, Nucleic Acids Res, 23:4407-4414 (1995)]. In addition, a PCR 
protocol, characterized by a highly shortened elongation time, can be used such that 
the amplification of short firagments is strongly favored thereby further reducing the 
25 numbo: and ttie average size of die amplicons. During the selective co-amplification 
of genomic firagments or in a subsequent PCR stq>, a first primer can be used that 
attaches a fiill promoter sequence (e.g., one deriving firom bacteriophage T7, T3 or 
SP6; supra) to the amplicons. The second strand may be synthesized by extension of 
a primer that contains a ribonucleotide residue at, for example, the penultimate 
30 position. Following PCR amplification, the primer sequences can be removed firam 



wo 00/66771 



PCT/EPOO/03904 



-52- 

this second strand by RNase digestion, and the resultant truncated strand transcribed 
with the aid of the first prima-. This procedure minimizes the common sequences that 
are connected to the target restriction fragments. 

5 EXAMPLE6 

cDNA Library Analyses - Transcription ProfiBng 

Diagnostic sequencing will, generally, be performed on a defined 
nucleic acid, one will know to what reference sequence the target nucleic acid 
corresponds. However, the re-sequencing methods according to the present invention 
10 can also be used to identify or classify certain sequences. In such experiments, the 

interrogated nucleic acid (e.g, a random clone of DNA) will typically correspond to an 
unknown portion of a (much) larger sample sequence or represent one out of a 
plurality of nucleic acids present in a biological san^>le, or a combination of both. The 
mass spectra derived fix)m the unknown nuclric add are compared to those known or 
1 S predicted for the related reference sequence(s), or portions thereof. Note that, in this 
type of experiments, some of the interrogated target sequences need not necessarily 
have their counteiparts in the reference sequences, and vice versa. It will be realized 
that sequence identification according to the present mefliods may, at the same time, 
reveal possible sequence variations. Litenogated sequences may thus be classified as 
20 identical to one of the database sequences, as a variant of such as a reference sequence 
or as novel in case no matching sequence is found 

It should be recognized that analyses that involve at least the four 
complemmtary mono-nucleotide specific cleavage reactions identify unknown 
sequences with a resolution essentially equal to sequence determination. At the same 
25 time, tiie MS-based methods described herein allow fast data acquisition and are 
amenable to hig^-throug^ut. Therefore, the present methods are usefiil to identify 
and catalogue nucleic acids at an unprecedented scale and speed. One ^plication 
consists of the analysis of cDNA libraries for the purpose of: (i) the assembly of 
unigene libraries (i.e the identification/removal of rq>Ucate clones), (ii) the 
30 identification of novel genes or novel variants of previously identified genes, and (iii) 
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transcription profiling. The speed and throughput of ttie present method should permit 
the processing of more clones and, hrace, a more in depth analysis of a cDNA hbrary. 

A variety of methods are known in the art for transcription profiling, 
i.e, the analysis of the transcription in both qualitative and quantitative terms. In one 

5 method, the expressed-sequence-tag (EST) approach, the mRNA population is 

assessed by partial sequencing of randomly selected cDNA clones. Global changes in 
gene-egression patterns are deduced fiom the EST ratios among two compared 
cDNA Ubraries [Lee N. et aL, Proc. Natl Acad, Set USA 92: 8303-8307 (1995)]. The 
methods described herein may be used to catalogue expressed genes with a similar 

10 level of resolution but considerably higher speed and throughput First, a library of 
unidiiectionally cloned cDNAs is constructed in a vector that permits transcription of 
the inserted sequences. Preferably, the 3*-aid of the cDNAs is located adjacent to the 
promoter. Tonplate for transcription can be prepared by amplification of the 
promoter-cDNA cassette using a pair of vector-^ecific primers. Alternatively, vector 

1 5 DNA is prq[>ared and cleaved at a restriction site within the vector close to the S *-end 
of the inserted cDNA (eg. -25 base-pairs). Preferably, the restriction site at which the 
templates are cleaved should have a low occurrence frequency within die cDNAs 
under study. Run-off transcripts, synthesized firom PCR product or digested vector, 
are diaracterized by a common 3'-end, consisting of vector sequences, which allows 

20 the isolation of fiill-length transcripts as described in Example 2. An alternative 

strategy involves treatment of the vector DNAs with a restriction reagent such that not 
only all templates are digested at the cDNA 5'-end but that a vast majority is also 
cleaved within the cDNA at some distance from the 3*-end {e.g. a few hundred base- 
pairs). The restriction reagent may be a single enzyme or a combination of two or 

25 more restriction enzymes. Ligation of an ad£q>tor to the digestion product(s) [see Vos 
P. et al. Nucleic Acids Res. 23: 4407-4414 (1995)] can be considered so as to obtain 
full-length transcripts with a common 3 '-end enabling their isolation as described in 
Example 2. However, transcripts that incorporate a biotin group at the 5 '-end may 
also be prepared [Hahner S. et al. Nucleic Acids Res,, 25: 1957-1964 (1997)], 

30 providing an alternative means for their immobilization. Digestion within the cDNAs 
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is an attractive option in that different partial cDNAs deriving bom the same 
transcript are made congruent by this procedure and thereby facile to identify. The 
full-length run-off transoipts are finally subjected to complementary sequence- 
specific cleavage reactions, and flie resultant digestion products analyzed by MS as 

5 disclosed herein. 

Those of sidll in the art will recognize the advantages of the transcript 
profiling metihod outlined above. Comparable to the EST approach, cDNAs are 
identified at the sequence-level, the ultimate level of resolution. Thus, while the 
method involves firagmentation of the interrogated nucleic acid, its level of resolution 

10 far exceeds that attained by fingerprinting techniques [Prashar Y. and Weissman S., 
Proc. Natl. Acad, Sci, USA 93: 659-663 (1996); Bachem C. etal, TTie Plant Journal 
9: 745-753 (1996); Ivanova N. and Belyavsky A., Nucleic Adds Res. 23: 2954-2958 
(1995); Liang P. and Pardee A., Science 257: 967-971 (1992)]. In contrast to 
hybridization-based approaches [Schena M. et al. Science 270: 467-470 (1995); 

15 Wodicka L. ei aL, Nature Biotechnology 15: 1359-1367 (1997)] the method can 

identify both known and previously unknown sequences. Also» it should prove faster 
then methods requiring gel-electrophoretic firactionation. 

EXAMPLE 7 

20 Whole-Genome Re-Sequencing 

In the past couple of years the technology for sequencing entire 
genomes, especially those of microorganisms, has come to maturity. More than 50 
microbial genomes are scheduled to be completed by the year 2000, and the benefits 
emerging fix>m this vast body of knowledge are nqiidly becoming clear [Clayton R. et 

25 al, Curr. Opinion Microbiol 1: 562-566 (1998)]. It seems clear that sequencing 
entire microbial genomes is becoming routine and that microbial genetics is entering 
the era of 'comparative genomics'. Knowledge of the complete genome sequence is 
the ultimate tool in phylogenetic analyses, allows gene/functional diversity studies, 
and fimdamentally changes the manner in which research is conducted in an 

30 organism. At the present time, a substantial portion of each new genome sequence has 
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no database match. One may expect to see a greater proportion of orthologous genes 
in the future, when the microbial species diversity is better represented. At that point, 
when most of the sequences generated will be similar to already known sequences, 
global genome analyses could be performed r^idly, accurately, and cost-efiectively 
5 using a re-sequencing strategy as desoibed herein rather than by de novo sequencing 
methods. Similar evolutions may be anticipated outside die bacterial genetics field 
where gmome projects for many (model) organisms are ongoing or have already been 
finished Drosophila melanogaster^ Caenorhabditis elegans^ human, mouse, 
Arabidopsis thaliana^ and rice). 
10 The methods of the present invention may be readily adapted to the re- 

sequencing of entire (bacterial) genomes or megabase nucleic acid regions. This may 
be accomplished with the use of a shotgun a^iproach that involves die sequence 
analysis of unselected subclones diat harbor random firagmoits according to the 
methods of die present invention. The assembly of all the independent, random 
1 5 sequences is fundamentally different fiom that in a ^fe now sequencing project 

ppleischmann R. et aL, Science 269: 496-512 (1995)] because of the availability of a 
reference sequence that serves as a scafifold The assembly into a single complete 
sequence comes down to matching each set of experimentally obtained spectra with a 
portion of the reference sequence. The computational approaches required to 
20 accompUsh this are similar to those that are needed for the analysis of cDNA libraries, 
outlined in Example 6. In both cases one does not know in advance the reference 
sequence, if at all existing, for a given interrogated target region. It should be noted, 
however, that the present shotgun qsproach might be even more demanding in tmns 
of computational power because of the undefined ends of the segments. At the same 
25 time, the algorithms must be csqiable of ms^ping the variations diat occur betwem the 
target and the reference sequence. It is expected that a shotgun abroach with its built- 
in redundancy (i.e., most sequences will be covered several-fold) should prove useful 
for the comprehensive comparison of a pair of related genomes. An altemative for the 
shotgun approach strategy consists of the analysis of clones from one or more libraries 
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of restriction enzyme fragments or the analysis of defined amplicons generated with 
locus specific primer pairs. 

While tibe present invention has been described in terms of the 
preferred embodiments, it is understood that variations and modifications will occur 
S to those skilled in the art Therefore, it is intended fliat Ae appended claims cover all 
such equivalent variations which come widiin the scope of the invention as claimed. 
All of the references dted herein are expressly incorporated by reference. 
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Table I: Detection of the twelve possible point mutations that can occur in DNA 
by the methods of the present invention. Each substitution is associated 
with the loss (- sign) and gain (+ sign) of a cleavage site. In addition, each 
mutation affects the mass of two digestion products as indicated. Mass 
5 difiTerences shown in bold face result firom the incoiporation of m^U in 

both transcripts (see text for details). 



Mutation 


RNaseTl 


RNaseU2 


(+) strand 1 (-)stxand 


(+)tiansciq>t | (-)tiansGr^ 


(+) transcript \ (-) transcript 


transitions 


A->G 


T->C 


+ 


-IDa 
-15 Da 




-1 Da 
-15 Da 


G->A 


C»T 




+lDa 
+15 Da 


+ 


+lDa 
+15 Da 


T->C 


A->G 


-1 Da 
*15Da 


+ 


-IDa 
-15 Da 




C->T 


a.>A 


+1 Da 
+15 Da 




+lDa 
+15 Da 


+ 


transversio 


ns 




T->G 


-24 Da 


+ 




+39 Da 
+25 Da 


C->A 


G->T 


+24 Da 




+ 


-39 Da 
-25 Da 


T->G 


A->C 


+ 


-24 Da 


+39 Da 
+25 Da 




G->T 


C->A 




+24 Da 


-39 Da 
-25 Da 


+ 


T->A 


A->T 


+23 Da 
49 Da 


-23 Da 
-9 Da 


+ 




A->T 


T->A 


-23 Da 
-9 Da 


+23 Da 
49 Da 




+ 


C->G 


G->C 


+ 




+40 Da 


^ODa 


G->C 


C->G 




+ 


-40 Da 


+40 Da 
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Table II: RNAse-A digestion products predicted for four different pGEM3-Zf(+) 
derived transcripts. The >3-mer fragments are ranked according to their 
molecular masses. The regular transcript was prepared with rNTP 
substrates. Transcripts that incoiporate dTMP, dUMP, or dCMP are 
denoted as dT-, dU-, or dC-transcript Fragments containing a S'- 
triphosphate (S'ppp-) are indicated. 







dtMrsttscHpl 


dC-tfsnse 


ripe 


Ffsgincnts 










expected mass 


ftigmetits 


expected hbss 
















CAT 


959.6 


TGC 


973.6 


TGC 


9593 


CCT 


903.S 


AAT 


983,6 


GAC 


9»3 


GAC 


9983 


CAT 


9433 


AGC 


998,6 


ATOC 


1302.8 


ATGC 


1288.8 


CAT 


9433 


ACC 


998,6 


AAGC 


1327.8 


AAGC 


13273 


AAT 


9833 


GAC 


998.6 


GAGC 


1343.8 


GAGC 


13433 


ACT 


9993 


ACT 


999.6 


AGGC 


1343.8 


AGGC 


1343.8 


C3CT 


10153 


GGC 


1014.6 


5^>PI^<5CGC 


1599.7 


TTCOC 


1594.9 


AGCT 


128B3 


OCT 


1015.6 


TTOGC 


1613.0 


5)i|i|^<XSGC 


1599.7 


AOCT 


12883 


OCT 


1013.6 


ATAGC 


16323 


ATAGC 


16183 


GGCT 


13043 


AAAT 


U12J 


GCTAC 


16483 


GOTAC 


1634.0 


AAAT 


1312.8 


AAGC 


13273 


TQTTTC 


1886.2 


TGTTTC 


1830.1 


GACT 


13443 


GAAT 


1328.8 


GAATTC 


1936.2 


GAATTC 


1908.1 


CACCT 


1521.9 


GAGC 


1343.8 


CTAATC 


1936a 


CTAATC 


1908.1 


GGCCT 


I6S03 


AGGC 


13433 


ATGGTC 


1952JI 


ATGGTC 


1924.1 


AGACT 


1674.0 


GACT 


13443 


TAGAGTC 


2281.4 


TAGAGTC 


2253.4 


CCACCT 


1867.1 


Sjipi^OCCSC 


1599.7 


GGGOATC 


2338.4 


GGGGATC 


2324.4 


CGAGCT 


1923.2 


AGACT 




TAAATAGC 


25943 


TAAATAGC 


2566.6 


CCAAGCT 


2252.4 


GGGGAT 




TATACTGTC 


28893 


TATAGTOTC 


2833.7 


GCAGGCAT 


25973 






TTOACTATTC 


31943 


TTGAOTATTC 


3123.9 


5)ipp-GGGCGAAT 


2893.5 






(5EQ1DN0: 22) 




(SEQIDNO:22) 




ACCCGGGGAT 
(SEOIDNa.23) 


32323 
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Table HI: RNase-A digestion products predicted for the dU- and dC-transcripts of 
the (+) and (-) sbands of the RNase-Tl coding region. Only the >3-mers 
are shown. Cleavage of the dU-transcript is C-specific. Likewise, the T- 
reaction is performed on the dC-transcript Two fragments, shown in 
5 italics, assume the occurrence of 3'*exteDded transcripts (refer to Example 

3). 



strand / CnreactifB 


1 IM+H)^ 


+) strand / T-nactioi 


,1 |M4^H)^ 


TTC 


904^ i 


\CT 


943.6 


TAC 


943.6 < 


SAT 


999.6 


TAC 


943»6 


cccr 


1192,7 


TAC 


943.6 


SGAT 


1344,8 


AAC 


982.6 


CCAAT 


1562,0 


AAC 


982.6 


GGCCr 


1594.0 


GAC 


998,6 


SA(jCr 


1634.0 


TATC^m* 


1153.7 


GAAACr 


1947,2 


TATV-pS' 


1233.7 


AOGAAT 


1947.2 


TTAC 


1233.7 


ACGAAGGT 


2637,6 


AATTC 


1562.9 


ACAACAACT 


2838,8 


S'p|i|i<SGGC 


1599.7 


5'ppp<3GGCGACr 


2853,5 


AAATAC 


1931.2 


ACCCACA(>^T 


3746,4 


GAAGAC 


20Q2.2 


(SEQ ID NO: 25) 




TGTGAGC 


2269.4 


CACGAAGACXKn* 


3890.4 


GAATGGC 


2308.4 


(SEQ ID NO: 26) 




GGTGAAAC 


2637«6 






TGTTGGATC 


2849.7 






GAAGGTTTTGATTrC (SEQ ID NO: 24) 


4723.8 






(.) strand /C-mrttan 


1 [M+Hl- 


(-) strand / T-iwtifM 


1 1 IM+HJ* 


TTC 


904.5 


CGT 


959.6 


TTC 


904.5 


CGT 


959.6 


GTC 


959.6 


CGT 


959.6 


AAC 


982.6 


AGT 


999.6 


5'ppp-GGGC 


1599.7 


AGT 


999.6 


AAAAC 


1641,0 


CCGG-OH3' 


1207.8 


AGTTTC 


1869.1 


CCCGG-OH3' 


1553,0 


GAATTC 


1908.1 


CKjAT 


1344.8 


AGAGAAATC 


2950.8 


G<K5T 


1360,8 


GTGAAGTTTATATC (SEQ ID NO: 27) 


4417,7 


GAAGT 


1674.0 


GTAGTAGGGAGAGC C5EQ ID NO: 28) 


4637,8 


CACCGT 


1867.1 


GTAG-inXfrroTA'rntjIXJI'GGGTAAGAATTGGATC 11123.7 


AAGAAT 


1987.2 


(SEQ ID NO: 20) 




CAAAACCT 


2509.6 






CCAACAGT 


2525.6 






ypn^GGGCGAAT 


2893.5 






AG(3GAGAGCr 


3328.0 






(SEQ ID NO: 29) 








CACAGAGAAAT 


3569.2 






(SEQ ID NO: 30) 
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Table IV: Spectral changes associated with single and multiple mutations in the 
RNase-Tl coding region. 



REACTION 


Awnrrvn nFFRRFNrR FRAGMElCrS 1 


NOVEL yRAGM» 1 




SEQUENCE Ii 








(+)/C 1 

(+)/T i 

( 

{-)/C i 
i 


mutatkm #1 fA->T on C+) atend/T-> 

[TAC 

VCCCACACAAAT 
SEQIDN0:2S) 

nAOTlVnvrATTTGTGTGGC^AAGAATTGGATC 
SEQIDNO:2q> 

SGCT 
\AGAAT 


A an t'-i 8» 

1233,7 
3746,4 

11123.7 

1360.8 
1987.2 


oDtttiiiinic 
not observed 


1194,7 

325.2 
3417.1 

11162,7 I 
3352,1 


<3-mer 
lot observed 


(+)/T 

«/c 


AAC 
FAC 

ACAACAACr 

GTACTTGTTGTATTTGTGTCKKn'AAGAATTGGATC 
(SEQIDNO:20) 

ACST 


»ConM 
982.6 
943.6 

2^8 

11.123J 

999.6 


not imiQue 
not oDiQiie 

not obscnfcd 
not unique 


1.947,2 

2.894.8 

1.288.8 
9.813.9 

943.6 


not observed 


(+>/C 

(-)/C 
(-)/T 


nmlirtkm «1 fA.>Tl» m smnd / TO 
GAAGOTrriOATTTC 
(SEQIDN0:2^ 

ACGAAGGT 
TTC 

CAAAACCT 
T 


>Aon(-) 
4723.8 

2637.6 

904.S 

2S09.6 
325,2 


notuniqye 
<3-nier 


4684.8 

1618.0 
1015.6 

943.6 

2838.8 


not observed 




mntatkm 14 f A->T OO f+> strand / T-> A cm M 








(+)/C 

(+)/T 

H/C 
(-)/T 


GAAGiil 1 flGATTTC 
(SEQ ID NO: 24) 

ACGAAGGT 

TTC 
T 


4723.8 

2637.6 

904.5 

325.2 
959.6 


not unique 

<3*nier 
not unique 


4684.8 

1288.8 
13443 

943.6 

1288.8 


notmuque 




mmmiai iS rAC->CG iiBmgtrnd/gr.>CGon(-)sli«nd! 






(+)/C 

(+)/T 
H/C 

(-)^T 


AAC 
TAC 

ACAACAACr 

GTAGL'Wi'roTA'iTitfTGTGGGTAAGAATTGGATC 
(SEQ ID NO: 20) 

ACT 
T 


982.6 
943.6 

2.838,8 

11.123.1 

999.6 
325.2 


DOS unique 
not unique 

r 

not unique 
<3-nier 


653.4 
1.288.8 

2.854,8 

1.288.8 
9.868.9 

1.288.8 


<3-mer 
not resolved 
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Table IV (continued) 



REACTION 


AFFECTED REFERENCE FRAGMENTS 


NOVEL FRAGM. 




SEQUENCE llM+HI^I Coomicnta 


1 Comncnls 




ButatioD 16 (AAA->CAG on (+) stnnd / XT 


r->CTG cm (•) strand] 






(+)/C 


AAATAC 


1931,2 




324,2 
1618.0 


<3-mer 


(+)/T 


ACCCACACAAAT 
(SEQIDN0:2^ 


3746.4 




3722.3 




(-)/c 


GTAbTlt?llUl'AlTX1j-imXKKn'AAGAATTGGATC 
(SEQID NO: 20) 


11123,7 




4104,4 
7108.3 




(-)/T 


AT 
T 

I 
GT 


654.4 
325,2 
325,2 
670.4 


<3-mei 
<3-iDer 
<3>iiier 
<3-incr 


943,6 
1015,6 




umtatioii #7 |AAT*>GCG od (+ ) strand / ATT->CGC oo (-) strand) 




(+)/C 


AATTC 


1S62.9 




669.4 
959,6 


<3-iner 


(+)/T 


CCAAT 
T 


1562.0 
325.2 


<3-mer 


1883.1 




(-)/c 


GTAGTTGTnn-ATTTGTGTGGGTAAGAATTGGATC 
(SEQID NO: 20) 


11123.7 




8904,3 
669,4 
1634.0 


<3-nier 


(-)/T 


AAGAAT 
1 

GGAT 


1987,2 
325,2 
1344,8 


<3-iner 


3601,2 






«8 (AACAAC, TAG -> AAG»AAT,TrC on (+) strand 


/ GTA»GTr,GTr -> CAA^Tr,CTr on (-) 


(+)/c 


AAC 
AA£ 
TAG 


982,6 
982,6 
943.6 


notuniqite 


2856,7 




(+)/T 


ACAACAACT 
ACGAAGGT 


2838.8 

2637.6 




2605.6 
325.2 
325.2 
2308.4 


<3-iner 
<3-iner 


(-)/C 


GTA(Jllt»n'rGTATilXyiUlWGTAAGAATTGGATC 
^EQ ID NO: 20) 


U123.7 




2237,4 
8888,3 




O/T 


CGT 
ACT 
GT 


III 


not unique 
notUDxpie 
<3-nier 


1947.2 
614.4 


<3-iner 
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Table V: U-specific cleavage of a 972 nucleotides long T7 transcript. The predicted 
digestion products, 222 in total, are grouped according to their 
composition. An asterisk indicates those peaks for which a companion 
<^clic phoq)hate reaction intermediate is observed (Figure 9). The largest 
S fragment is absent from the obtained spectrum; a few other cleavage 

products appear as minor peaks and are labeled *weak\ 



C!oninnsition 


M+H 


Length 


Number 


Remarks 


T 


325^ 


1 


47 




CT 


614,4 


2 


11 




AT 


6S4»4 


2 


14 




GT 


670,4 


2 


15 






903,5 


3 


4 




ACT 


943.6 


3 


3 




CGT 


959,6 


3 


7 




AjT 


983,6 


3 


5 




ACT 


999,6 


3 


1 




GjT 


1015,6 


3 


4 




CjT 


1192,7 


4 


2 




ACjT 


1232,7 


4 


5 




CjCT 


1248,7 


4 


4 




AjCT 


1272,8 


4 


3 




ACX5T 


1288,8 


4 


6 




CGjT 


1304,8 


4 


5 




A,T 


1312,8 


4 


1 


weak 


AjCT 


1328,8 


4 


1 


weak 


AGjT 


1344,8 


4 


5 




ACjT 


1521,9 


5 


1 




C3GT 


1537.9 


5 


2 




AAT 


1562,0 


5 


2 




AC^CT 


1578,0 


5 


2 




C,G,T 


1594,0 


5 


7 


* 


AjCCT 


1618,0 


5 


1 


weak 


ACXj^T 


1634,0 


5 


3 




CG3T 


1650,0 


5 


3 




A3GT 


1658,6 


5 


2 




AAT 


1674,0 


5 


2 




G4T 


1706,0 


5 


1 




C5T 


1771,1 


6 


1 




C4CT 


1827,1 


6 


1 




AC5GT 


1867,1 


6 


2 


* 


C,G,T 


1883,1 


6 


2 




A,C,T 


1891,2 


6 


1 
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TableV (continued) 



AC fm-T 




6 


2 


* 






5 










g 










/ 






C4O2T 




7 










7 






A 0 rf7» T 
AC5CJ2T 




•7 










/ 






A p«i 
A2C2G2T 


2252,4 


/ 






^ ^ t* 
C2G4T 


22o4,4 


/ 






AjCGaT 


2292,4 


/ 






AjCGjT 


2308,4 


7 






AOG4T 


2324,4 


7 






ACsOT 


2445,5 


8 


J 




AAG3T 


2597,6 


8 






A4CG2T 


2621,6 


8 






A2CG4T 


2633,6 


0 
0 


■ 




A4C3GT 


2854,8 


9 








2878,8 


A 

9 






AjCjGjT 


2886,8 


9 






A2CXj4T(5*iipp-) 


2893,6 


0 
0 


j 




A3C2G3T 


2926,8 


A 

9 


J 




C,GT 


2983,9 


lU 




WCSK 


A2C3G2T 


3119,9 


1 A 

lU 






A ^ ^ MN 


321 0,0 


in 










111 






A 0 n T 


•^979 A 


in 






A n T 


0900 A 


in 
ill 






A5C5T 


liin 1 
J4i /,l 


1 1 






AjCjGiT 


3529,2 


11 






A^CXJjT 


3625,2 


11 






ACjG^T 


3938,4 


12 






A^C^jT 


4043,5 


13 






AjCjGjT 


4139,6 


12 






AAG4T 


4219,6 


13 






A4C2G5T 


4291,6 


13 






A,C,G,T 


4661,9 


15 






A5C4G5T 


4854,0 


15 






A,C,G4T 


5536,4 


17 






A,QG,T 


6106,8 


19 






A4C^,3T 


8154,0 


25 






AiaQGioT 


103704 


32 




not observed 






2=222 
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We Claim : 

1 . A method for sequence analysis of one or more target nucleic acids 
for which a refermce nucleic add sequence is known, ttie method comprising the 

S steps of: 

(a) deriving from one or more biological samples the one or 
more target nucleic acids; 

(b) subjecting the one or more target nucleic acid to 
complementary cleavage reactions with one or more cleavage reagents thereby 

10 generating cleavage products; 

(c) performing mass spectioscopical analysis on the cleavage 
products generated in step (b) to obtain one or more mass spectra; and 

(d) comparing the one or more mass spectra of the cleavage 
products of the one or more target nucleic acids obtained in step c) with a known or 

1 S predicted mass spectra for the reference nucleic acid sequence, and deducing 

therefrom, by systematic computational analysis, all or part of the nucleotide sequence 
of the one or more target nucleic acids, and comparing the deduced nucleic acid 
sequence with the reference nucleic acid to determine whether the one or more target 
nucleic acids have fhe same sequence or a different sequence from the reference 

20 nucleic acid. 

2. The method of claim 1 wherein the nucleic acid sequence 
difference that is detemiined in step d) is a deletion, substitution, insertion or 
combinations thereof. 

25 

3. The method of claim 1 wherein the biological sample is derived 
from organism selected from the group consisting of eukaryotes, prokaryotes, and 
viruses. 



30 



4. 



The method of claim 1 wherein the target nucleic acid is selected 
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from the ffoup consisting of single stranded DNA, double stranded DNA, cDNA» 
single stranded RN A, double stranded RNA, DN A/RN A hybrid, and DN A/RNA 
mosaic nucleic acid. 

S 5. The method ofclaim 4 wherein the one or more target nucleic 

acids is derived by one or more consecutive anq[>lification procedures selected fiom 
the group consisting of in vivo cloning, polymerase chain reaction (PCR), reverse 
transcription followed by the polymerase chain reaction (RT-PCR), strand 
displacement amplification (SD A), and transcription based processes. 

10 

6. The method of claim 5 wherein the one or more ampUfied target 
nucleic acids are RNA transcripts generated fiom a single stranded or a double 
stranded target nucleic acid by a process oonqirising the stq>s of: 

(a) linking operatively an expression control sequences to the 
1 S one or more target nucleic acids; and 

(b) transcribing one or both strands of the one or more target 
nucleic acid of step a) using one or more RNA polymerases that recognize the 
transcription control sequence on the one or more target nucleic acids. 

20 7. The method of claim 6 wherein the transcription control sequence 

is selected from the group consisting of an eukaryotic transcription control sequence, 
a prokaryotic transcription control sequence, and a viral transcription control 
sequence. 

25 8. The method of claim 7 wherein the prokaryotic transcription 

control sequence is selected from the group consisting of T3, T7, and SP6 promoters. 

9. The method of claim 8 wherein the RNA polymerases which 
utilize the T3, T7, or SP6 promoters are either wild type or mutant RNA polymerases. 
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the mutant polymerases being capable of incorporating into the RNA transcript non- 
canonical substrates with a 2'-substituent other than a hydioxyl group. 

10. The method of claims 4, S, 6, 7, 8 or 9 wha:ein the target nucleic 

S acid is dmved using one or more modified nucleoside triphosphates selected fiom the 
group consisting of mass modified deoxynucleoside triphosphates, mass modified 
dideoxynucleoside triphosphates, and mass modified ribonucleoside triphosphates. 

1 1 . The method of claim 10 wherein the modified nucleoside 

10 triphosphates are modified on the base, sugar, and/or die phosphate moiety, and are 
introduced through an enzymatic stq>, chemically, or a combination of both. 

12. The m^hod of claim 10 whoein the modification consists of 2'- 
substituents other dum a hydroxyl group on the nucleotide triphosphates. 

15 

13. The method of claim 10 wherein the modification consists of 
phosphorothioate intemucleoside linkages or phosphoroduoate intemucleoside 
linkages further reacted with an alkylating reagent 

20 14. The method of claun 10 wherdn the modification consists of a 

methyl group on C5 of the uridine-5 '-monophosphate subunits. 

15. The method of claim 10 wherein the modifications alter the 
specificity of cleavage by the one or more cleavage reagents and/or the mass and/or 

25 the length of the cleavage products. 

16. The method of claims 11, 12, 13 or 14 wherein the modifications 
alter the specificity of cleavage by die one or more cleavage reagents and/or the mass 
and/or the lengOi of the cleavage products. 
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17. The method of claim 1 wherein the complementary cleavage 
reactions are selected from the group consisting of enzymatic cleavage, chemical 
cleavage, and physical cleavage. 

5 

18. The method of claim 17 wherein the complemmtary cleavage 
reactions are characterized by a relaxed mono-nucleotide, mono-nucleotide, relaxed 
di-nucleotide, or di-nucleotide specificity. 

10 19. The method of claim 1 7 or claim 1 8 wherein the chemical 

cleavage reaction consists of an alkaU treatment, and the one or more target nucleic 
acids are RNA transcripts that have incorporated non-canonical substrates with a 2 
substituent oHier than a hydroxyl groiqp. 

IS 20. The method of claim 17 or claim 18 wherein the one or more 

target nucleic acids are subjected to enzymatic cleavage reaction using one or more 
enzymes selected fit>m the groiq> consisting of endonucleases and exonucleases. 

21 . The method of claim 20 wherein the one or more target nucleic 
20 acids are subjected to enzymatic cleavage reaction using one or more endonucleases, 

selected from the group consisting of restriction enzymes, RNA endonucleases, DNA 
endonucleases and non-specific phosphodiesterases. 

22. The method of claim 2 1 wherein the one or more endonucleases 
25 are one or more selective or non-selective KNA endonucleases and the one or more 

target nucleic acids are transcripts. 

23. The method of claun 22 wherein the one or more selective or non- 
selective RNA endonuclease is selected from the group consisting of the G-specific Tj 
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ribonuclease, the A-spedfic ribonuclease, the A/U specific phyM ribonuclease, the 
U/C specific ribonuclease A, the C-q)ecific chicken liver ribonuclease (RNaseCLS), 
and cusativin. 

S 24. The method of claim IS wherein the one or more target nucleic 

adds are phosphorothioate-modified single-stranded DNA or RNA and &e 
endonuclease is nuclease PI. 

25. The method of claim 1 whorein the mass spectroscopical analysis 
10 of the nucleic acid fifagments is performed using a mass spectroscopic method 

selected fiom the group consisting of Matrix-Assisted Laser Desorption/Ionization- 
Time-of-flight (MALDI-TOF), Electiospiay-Ionization (ESQ, and Fourier Transfomi- 
lon CyclotEon Resonance (FT-ICR). 

IS 26. A method for scoring known nucleotide sequence variations of one 

or more target nucleic acids for which a known reference nucleic acid sequence is 
available, the method comprising the stq)s of: 

(a) ' deriving fiom one or more biological samples the one or 

more target nucleic acids; 
20 (b) subjecting the one or more target nucleic acids to 

complementary cleavage reactions with one or more cleavage reagents thereby 

generating cleavage products; 

(c) performing mass spectroscopical analysis on the cleavage 

products obtained in stq> (b) to obtain two or more mass spectra; and 
2S (d) comparing die one or more mass spectra of the cleavage 

products of the one or more target nucldc adds obtained in step c) witii the known or 

predicted mass spectra for the ref^ence nucldc add sequence, and scoring therefiom, 

by systonatic computational analysis the presence or absence of the known nucleotide 

sequence variations of flie one or more target nucldc acids. 
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27. The method of claim 26 wherein tiie nucleic acid sequence 
variation of the target nucleic acid scored in step d) is a deletion, a substitution, 
insertion or combinations theteof 

5 

28. The method of claim 26 wherein flie one or more target nucleic 
adds are derived firom biological samples selected fiom tiie groiq> consisting of an 
eukaiyote, a prokaryote, and a virus. 

1 0 29. The method of claim 26 or claim 27 wherein the one or more 

target nucleic acids are selected from the group consisting of a single stranded DNA» a 
double stranded DNA, a single stranded RNA» a double stranded RNA, and a 
DNA/RNA hybrid. 

IS 30. The method of claim 29 wherein the one or more target nucleic 

acids are derived by one or more consecutive amplification procedures selected from 
the group consisting of m vivo cloning, the polymerase chain reaction (PCR), reverse 
transcription foUowed by the polymoase chain reaction (RT-PCR), strand 
displacement amplification (SDA), and transcription based processes. 

20 

3 1 - The method of claim 29 wherein the one or more amplified target 
nucleic acids are RNA transcripts generated by a process comprising the steps of: 

(a) linking operatively a transcription control sequences to the 
one or more target nucleic acids; and 
25 (b) transcribing one or both strands ofthe one or more target 

nucleic acids using one or more RNA polymerases that recognize the transcription 
control sequence on the one or more target nucleic adds. 



32. The method of claim 31 wherein the transcription control sequence 
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ls selected from the group consisting of an eukaiyotic transcription control sequence, 
a prokaryotic transcription control sequence, and a viral transcription control 
sequence. 

S 33. The inetfaodofclaim 32 wherein the prokaryotic transcription 

control sequence is selected fiom the group consisting of T3» T7, or SP6 promoters. 

34. The method of claim 33 wherein the RNA polymerases which 
utilize the T3, T7 or SP6 promoters are either wild type or in a mutant form capable of 

1 0 incorporating non-canonical substrates with a 2'-substituent other than a hydroxyl 
groiqp. 

35. The mefliod of claim 30, 31 or 32 wherein the one or more target 
nucleic acids are derived using a modified nucleoside triphosphates selected fiom the 

1 5 gjxnxp ccmsisting of mass modified deoxynucleoside triphosphates, mass modified 
dideoxynucleoside triphosphates, and mass modified ribonucleoside riphosphates. 

36. The method of claim 35 whmin the modified nucleoside 
triphosphates are modified on the base, sugar, and/or the phosphate moiety, and are 

20 introduced through an enzymatic step, chemically, or a combination of both. 

37. The mefliod of claim 35 wherein the modification consists of 2^* 
substituents other than a hydroxyl gfovp on the nucleoside triphosphates. 

38. The method of claim 35 who^ the modification consists of 
25 phosphorothioate intemucleoside linkages or phosphorothioate intemucleoside 

linkages Anther reacted with an alkylating reagent 

39. The method of claim 35 wherein the modification consists of a 
methyl group on C5 of the uridine-5 '-monophosphate subunits. 
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40. The mefliod of claim 35 wherein the modificalions alter the 
specificity of cleavage by the one or more cleavage reagents and/or the mass and/or 
the length of the cleavage products. 

5 

41. The method of claim 36, 37, 38, 39 or 40 wherein the 
modifications alter the specificity of cleavage by the one or more cleavage reagents 
and/or the mass and/or the length of the cleavage products. 

10 42. The method of claim 26 wherein the complementary cleavage 

reactions are firom the groiq> consisting of enzymatic cleavage, chemical cleavage, and 
physical cleavage. 

43. The mediod of claim 42 wherein the complonentary cleavage 

1 5 reactions are characterized by a relaxed mono-nucleotide, mono-nucleotide, relaxed 
di-nucleotide, or di-nucleotide specificity. 

44. The method of claim 42 or claim 43 wherein the chemical 
cleavage consists of an alkali treatment, and the one or more target nucleic acids are 

20 RNA transcripts that incorporate non-canonical substrates with a 2'-substituent other 
than a hydroxyl group. 

45. The method of claim 42 or claim 43 wherein the one or more 
target nucleic acids are subjected to enzymatic cleavage reaction using one or more 

25 enzymes selected fiom the group consisting of endonucleases and exonucleases. 

46. The method of claim 45 wherein the one or more target nucleic 
acids are subjected to enzymatic cleavage reaction using one or more endonucleases, 
selected fipom the group consisting of restriction enzymes, RNA endonucleases, DNA 
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endonucleases, and non-specific phosphodiesterases. 

47. The method of claim 46 wherein the one or more endonucleases 
are the one or more selective or non-selective RNA endonucleases and one or more 

S target nucleic acids are transcripts. 

48. The method of claim 46 or claim 47 wherein the one or more RNA 
endonucleases are selected fix)m the group consisting of the G-specific T, 
ribonuclease, the A-specific rihonuclease, the A/U specific phyM ribonuclease, the 

10 U/C specific ribonuclease A, the C-specific chicken liver ribonuclease (RNaseCL3)» 
and cusativin. 

49. The mefliod of claim or claim wherein the one or more target 
nucleic adds are phosphorothioate- modified single-stranded DNA or RNA and the 

IS endonuclease is nuclease PI . 

50. The method of claim 26 wherein the mass spectroscopical analysis 
of the nucleic acid fragments is performed using a mass spectroscopic method 
selected firom the groiq> consisting of Matrix-Assisted Laser Desorption/Ionization- 

20 Time-of-flight (MALDI-TOF), Electro^ray-Iomzation (ESI), and Fourier Transform- 
Ion Cyclotron Resonance (FT-ICR). 

51. A method for detomining the sequence of one or more target 
nucleic acids, the method con^irising the stq>s of: 

2S (a) dmving from a biological sample material the one or more 

target nucleic acids; 

(b) subjecting the one or more target nucleic acids to 
complementary cleavage reactions with one or more cleavage reagents thereby 
generating cleavage products; 
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(c) perfonning mass spectroscopical analysis on the cleavage 
pioducts obtained in stq> (b) to obtain a mass spectra of the nucleic acid fragments; 
and 

(d) deducing therefrom the nucleotide sequence of the one or 

S more target nucleic acids by systematic computational analysis on die mass spectra of 
the cleavage products obtained in step (c). 

52. A method for genome wide genotyping of one or more known or 
unknown target nucleic acids, the method comprising the steps of: 

10 (a) deriving from one or more biological samples the one or 

more target nucleic acids; 

(b) subjecting the one or more target nucleic acids to one or 

more specific cleavage reactions with one or more cleavage reagents thereby 

generating cleavage products; 
1 5 (c) performing mass spectroscopical analysis on the cleavage 

products obtained in st^ (b) to obtain one or more mass spectra; and 

(d) comparing the mass spectra of the cleavage products 

obtained in step c) with reference mass spectra, and diagnosing therefrom genetically 

relevant nucleic acid sequence variations of the one or more known or unknown target 
20 nucleic acids. 

53. A method for identifying one or more target nucleic acids in one or 
more biological samples, the method comprising the steps of: 

(a) deriving from one or more biological samples the one or 
25 more target nucleic adds; 

(b) subjecting the one or more target nucleic acids of step (a) 
to one or more specific cleavage reactions with one or more cleavage reagents thereby 
generating cleavage products; 

(c) performing mass spectroscopical analysis on the cleavage 
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pioducts obtained in stq> (d) to obtain one or more mass spectra; and 

(d) coiqpaiing the one or more mass sp ectra obtained in step c) 
with each other or witib a plurality of mass spectra of reference nucleic acids, and 
deducing therefrom the identity of the one or more target nucleic adds. 

5 

54. The method of claim 53 wherein the one or more target nucleic 
adds are cDNA. 

55. The method of claim 54 wherein the method is used to determine 
10 the expression profile in a biological sample. 

56. A kit for sequence analysis of one or more target nucleic acids 
using mass spectroscopy, the kit comprising: 

(a) one or more sets of reference nucleic acids for which the 
1 5 nucldc add sequence is knowi^ 

(b) one or more nucleotide triphosphates; 

(c) one or more polymerases; 

(d) one or more nucldc acid cleaving agents; and 

(e) computer software for comparing the mass spectra of the 
20 one or more target nucldc acid with the mass spectra of the reference nucleic acid and 

deducing therefrom the nucleic acid sequmce of the target nucleic acid. 
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SEQUENCE LISTING 



<110> METHEXIS MV 

<120> SEQUENCING BY A COMBINATION OP MONONUCLEOTIDE-SPBCIFIC 
DIGESTI(»9 AND MASS SPBCTROHETRY 

<130>. 29314/35410A 

<140> 
<141> 

<150> 60/131,984 
<151> 1999^04-30 

<160> 30 

<170> Patentin Ver. 2,1 

<210> 1 

<211> 120 

<212> DNA 

<213> Homo sapiens 

<220> 

<223> exon 5 o£ human p53 
<400> 1 

tactcccctg ccctcaacaa gatgttttgc caactggcca agacctgccc tgtgcagctg 60 
tgggttgatt ccacaccccc gcccggcacc cgcgtccgcg ccatggccat ctacaagcag 120 

<210> 2 
<211> 40 
<212> DNA 

<213> Artificial Sequence 
<220> 

<223> Description of Artificial Sequence: synthetic 
<220> 

<223> pGEM3-Z£(^) derived nucleotide 
<400> 2 

gtaaaacgac ggccagtgaa ttgtaatacg actcactata 40 



<210> 3 
<211> 972 
<212> DNA 

<213> Artificial Sequence 
<220> 

<223> Description of Artificial Sequence: synthetic 
<220> 

<223> pGEM3-Zf ( + ) derived nucleotide 
<400> 3 

gggcgaattc gagctcggta cccggggatc ctctagagtc gacctgcagg catgcaagct 60 
tgagtattct atagtgtcac ctaaatagct tggcgtaatc atggtcatag ctgtttcctg 120 
tgtgaaattg ttatccgctc acaattccac acaacatacg agccggaagc ataaagtgta 180 
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aagcctgggg tgcctaatga gtgagctaac 
ctttccagtc gggaaacctg tcgtgccagc 
gaggcggttt gcgtattggg cgctcttccg 
tcgttcggct gcggcgagcg gtatcagctc 
aatcagggga taacgcagga aagaacatgt 
gtaaaaaggc cgcgttgctg gcgtttttcc 
aaaatcgacg ctcaagtcag aggtggcgaa 
ttccccctgg aagctccctc gtgcgctctc 
tgtccgcctt tctcccttcg ggaagcgtgg 
tcagttcggt gtaggtcgtt cgctccaagc 
ccgaccgctg cgccttatcc ggtaactatc 
tatcgccact ggcagcagcc actggtaaca 
ctacagagtt cttgaagtgg tggcctaact 
tctgcgctct gc 



-2- 

tcacattaat tgcgttgcgc tcactgcccg 240 
tgcattaatg aatcggccaa cgcgcgggga 300 
cttcctcgct cactgactcg ctgcgctcgg 360 
actcaaaggc ggtaatacgg ttatccacag 420 
gagcaaaagg ccagcaaaag gccaggaacc 480 
ataggctccg cccccctgac gagcatcaca 540 
acccgacagg actataaaga taccaggcgt 600 
ctgttccgac cctgccgctt accggatacc 660 
cgctttctca tagctcacgc tgtaggtatc 720 
tgggctgtgt gcacgaaccc cccgttcagc 780 
gtcttgagtc caacccggta agacacgact 840 
ggattagcag agcgaggtat gtaggcggtg 900 
acggctacac tagaagaaca gtatttggta 960 

972 



<210> 4 
<211> 131 
<212> DNA 

<213> Artificial Sequence 
<220> 

<223> Description o£ Artificial Sequence: synthetic 
<220> 

<223> PGR products and transcripts 
<400> 4 

taatacgact cactataggg cgacttcacg aagacggtga aactgttgga tccaattctt 60 
acccacacaa atacaacaac tacgaaggtt ttgatttctc tgtgagctct ccctactacg 120 
aatggcctat c 131 



<210> 5 
<211> 134 
<212> DNA 

<213> Artificial Sequence 
<220> 

<223> PGR products and transcripts 
<220> 

<223> Description of Artificial Sequence: synthetic 
<400> 5 

taatacgact cactataggg cgaattcgta gtagggagag ctcacagaga aatcaaaacc 60 
ttcgtagttg ttgtatttgt gtgggtaaga attggatcca acagtttcac cgtcttcgtg 120 
aagtttatat ccgg 134 



<210> 6 
<211> 45 
<212> DNA 

<213> Artificial Sequence 
<220> 

<223> Description of Artificial Sequence: synthetic 
<220> 

<223> reference nucleotide 



<400> 6 

ggatccaatt cttacccaca caaatacaac aactacgaag gtttt 



45 



wo 00/66771 



PCT/EPOO/03904 



-3- 



<210> 7 
<211> 45 
<212> DNA 

<213> Artificial Sequence 
<220> 

<223> Description of Artificial Sequence: synthetic 
<220> 

<223> mutant 1 
<4D0> 7 

ggatccaatt ctttcccaca caaatacaac aactacgaag gtttt 45 

<210> 8 
<211> 45 
<212> DNA 

<213> Artificial Sequence 
<220> 

<223> Description of Artificial Sequence: synthetic 
<220> 

<223> mutant 2 
<400> 8 

ggatccaatt cttacccaca caaatacaac aagtacgaag gtttt 45 

<210> 9 
<211> 45 
<212> DNA 

<213> Artificial Sequence 
<22a> 

<223> Description of Artificial Sequence: synthetic 
<220> 

<223> mutant 3 
<400> 9 

ggatccaatt cttacccaca caaatacaac aactacgatg gtttt 45 

<210> 10 
<211> 45 
<212> WA 

<213> Artificial Sequence 
<220> 

<223> Description of Artificial Sequence: synthetic 
<220> 

<223> mutant 4 
<400> 10 

ggatccaatt cttacccaca caaatacaac aactacgtag gtttt 45 



<210> 11 
<211> 45 
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<212> DNA 

<213> Artificial Sequence 
<220> 

<223> Description of Artificial Sequence: synthetic 
<220> 

<223> mutant 5 
<4Q0> 11 

ggatccaatt cttacccaca caaatacaiac acgtacgaag gtttt 45 



<210> 12 
<211> 45 
<212> DNA 

<213> Artificial Sequence 
<220> 

<223> Description of Artificial Sequence: synthetic 
<220> 

<223> mutant 6 
<400> 12 

ggatccaatt cttacccaca ccagtacaac aactacgaag gtttt 45 

<210> 13 
<211> 45 
<212> DNA 

<213> Artificial Sequence 
<220> 

<223> Description of Artificial Sequence: synthetic 

<220> 

<223> mutant 7 
<400> 13 

ggatccgcgt cttacccaca caaatacaac aactacgaag gtttt 45 

<210> 14 
<211> 45 
<212> DNA 

<213> Artificial Sequence 
<220> 

<223> Description of Artificial Sequence: synthetic 
<220> 

<223> mutant 8 
<400> 14 

ggatccaatt cttacccaca caaatacaag aatttcgaag gtttt 45 

<210> 15 
<211> 13 
<212> DNA 

<213> Artificial Sequence 
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<220> 

<223> Description of Artificial Sequence: synthetic 



<400> 15 
ctagcccccg ate 



13 



<210> 16 
<211> 25 
<212> DNA 

<213> Artificial Sequence 
<220> 

<223> Description of Artificial Sequence: primer 
<400> 16 

ccggatataa acttcacgaa gacgg 25 



<210> 17 
<211> 26 
<212> DNA 

<213> Artificial Sequence 
<220> 

<223> Description of Artificial Sequence: primer 



<210> 18 
<211> 37 
<212> DNA 

<213> Artificial Sequence 
<220> 

<223> Description of Artificial Sequence: primer 
<400> 18 

taatacgact cactataggg cgacttcacg aagacgg 37 



<210> 19 
<211> 41 
<212> DNA 

<213> Artificial Sequence 
<220> 

<223> Description of Artificial Sequence: primer 



<400> 17 

gataggccat tcgtagtagg gagagc 



26 



<400> 19 

taatacgact cactataggg cgaattcgta gtagggagag c 



41 



<210> 20 
<211> 35 
<212> DNA 



<213> Artificial Sequence 



<220> 

<223> Description of Artificial Sequence: synthetic 



wo 00/66771 



-6- 

<220> 

<223> reference fragment 
<400> 20 

gtagttgttg tatttgtgtg ggtaagaatt ggatc 



<210> 21 
<211> 15 
<212> PN?k 

<213> Artificial Sequence 
<220> 

<223> Description of Artificial Sequence: synthetic 
<220> 

<223> polymorphic 15-mer fragment 
<400> 21 

aaauceuiaac cuucg 



<210> 22 
<211> 10 
<212> ONA 

<213> Artificial Sequence 
<220> 

<223> RNAse-A digestion products 
<220> 

<223> Description of Artificial Sequence: synthetic 

<400> 22 
ttgagtattc 



<210> 23 
<211> 10 
<212> DMA 

<213> Artificial Sequence 
<220> 

<223> RNAse-A digestion products 
<220> 

<223> Description of Artificial Sequence: synthetic 

<400> 23 
acccggggat 



<210> 24 
<211> 15 
<212> DNA 

<213> Artificial Sequence 
<220> 

<223> RNAse-A digestion products 
<220> 

<223> Description of Artificial Sequence: synthetic 



wo 00/66771 



PCT/£P0Qm3904 



.7- 



<400> 24 

gaaggttttg atttc 



15 



<210> 25 
<211> 12 
<212> DNA 

<213> Artificial Sequence 
<220> 

<223> KN2^e-A digestion products 
<22D> 

<223> Description of Artificial Sequence: synthetic 
<400> 25 

acccacacaa at 12 



<210> 26 
<211> 12 
<212> DNA 

<213> Artificial Sequence 
<220> 

<223> RNAse-A digestion products 
<220> 

<223> Description of Artificial Sequence: synthetic 



<210> 27 
<211> 14 
<212> DNA 

<213> Artificial Sequence 
<220> 

<223> RMAse-A digestion products 

<220> 

<223> Description of Artificial Sequence: synthetic 
<400> 27 

gtgaagttta tatc 14 



<210> 28 
<211> 14 
<212> DNA 

<213> Artificial Sequence 
<220> 

<223> RNAse-A digestion products 
<22Q> 

<223> Description of Artificial Sequence: synthetic 



<400> 26 
cacgaagacg gt 



12 



<400> 28 
gtagtaggga gage 



14 



wo 00/66771 



PCT/EPOO/03904 



-8- 

<210> 29 
<211> 10 
<212> DNA 

<213> Artificial Sequence 
<220> 

<223> RNAse-A digestion products 
<220> 

<223> Description of Artificial Sequence: synthetic 
<400> 29 

agggagagct 10 



<210> 30 
<211> 11 
<212> DNA 

<213> Artificial Sequence 
<220> 

<223> RNAse-A digestion products 



<220> 

<223> Description of Artificial Sequence: synthetic 



<400> 30 
cacagagaaa t 
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